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Abstract 

We study a model of n one-dimensional non-intersecting Brownian motions with two 
prescribed starting points at time t — and two prescribed ending points at time t = 1 
in a critical regime where the paths fill two tangent ellipses in the time-space plane as 
n —¥ oo. The limiting mean density for the positions of the Brownian paths at the time of 
tangency consists of two touching semicircles, possibly of different sizes. We show that in 
an appropriate double scaling limit, there is a new familiy of limiting determinantal point 
processes with integrable correlation kernels that are expressed in terms of a new Riemann- 
Hilbert problem of size 4x4. We prove solvability of the Riemann-Hilbert problem and 
establish a remarkable connection with the Hastings-McLeod solution of the Painleve II 
equation. We show that this Painleve II transcendent also appears in the critical limits 
of the recurrence coefficients of the multiple Hermite polynomials that are associated with 
the non-intersecting Brownian motions. Universality suggests that the new limiting kernels 
apply to more general situations whenever a limiting mean density vanishes according to 
two touching square roots, which represents a new universality class. 

Keywords: non-intersecting Brownian motion, determinantal point process, correlation 
kernel, Riemann-Hilbert problem, Deift-Zhou steepest descent analysis, Painleve II equation, 
multiple Hermite polynomial. 
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1 Introduction 

In recent years the model of non-intersecting Brownian motions has been studied in various 
regimes, see e.g. El EH [23 123 Ell E3 HI], where many connections with determinantal point 
processes and random matrix theory were found, see also [32l |37l [36l 04] for non-intersecting 
Bessel paths and Brownian excursions. Their discrete counterparts, the non-intersecting random 
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Figure 1: Non-intersecting Brownian motions with two starting and two ending positions in case 
of (a) large, (b) small, and (c) critical separation between the endpoints. Here the horizontal axis 
denotes the time, t G [0, 1], and for each fixed t the positions of the n non-intersecting Brownian 
motions at time t are denoted on the vertical line through t. Note that for n — > oo the positions 
of the Brownian motions fill a prescribed region in the time-space plane, which is bounded by 
the boldface lines in the figures. Here we have chosen N = n = 20 and p\ = p 2 = 1/2 in each of 
the figures, and (a) <zi = — a 2 = 1, b\ = —b 2 = 0.7, (b) a\ = — a 2 = 0.4, b\ = — 6 2 = 0.3, and (c) 
a\ = —a-i = 1, b\ = —b 2 = 1/2, in the cases of large, small and critical separation, respectively. 

walks, have important connections with tiling and random growth models, see e.g. [51 IT21 1^51 [2^1 

021- 
Iri this paper we consider one-dimensional non-intersecting Brownian motions with two pre- 
scribed starting points at time t = and two prescribed ending points at time t = 1. We assume 
that the number of paths that leave from the topmost (bottommost) starting point is the same 
as the number of paths that arrive at the topmost (bottommost) ending point. As the number 
of paths increases and simultaneously the overall variance of the Brownian transition probability 
decreases, we may create various situations that are illustrated in Figure [T] 

In case of large separation of the starting and ending points we have a situation as in Fig- 
ure HJa). The paths are in two disjoint groups, where one group of paths goes from the topmost 
starting point to the topmost ending point, and the other group goes from the bottommost 
starting point to the bottommost ending point. In the large n limit the paths fill out two disjoint 
ellipses. 

In case of small separation of the starting and ending points we have a situation as in Fig- 
ure [TJb) . Here the two groups of paths that emanate from the two starting points merge at a 
certain time, stay together for a while, and separate at a later time. In the large n limit the paths 
fill out a region that is bounded by a more complicated curve with two cusp points. The critical 
behavior at the cusp point is known to be described by the Pearcey process, see [SJ EJ HI] ■ 
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In the transitional case of critical separation the paths fill out two ellipses that are tangent 
at a critical time as shown in Figure QJc). The case of critical separation was already considered 
by the first two authors [15] . but at the non-critical time. Here we consider the behavior at the 
critical time. Note that the tangent point of two ellipses is called a tacnode in the classification 
scheme of singular points of algebraic curves, whence the title of this paper; see also [TT] where 
another model is analyzed with a tacnode, but with markedly different properties. 

The phase transition at critical separation can already be observed at times different from 
the critical time t CT it. In [T!5] a connection was found with the Hastings-McLeod solution of the 
Painleve II equation 

q"( S ) = S q( S ) + 2q 3 ( S ). (1.1) 

Here the prime denotes the derivative with respect to s. The Hasting-McLeod solution [25] is 
the special solution q(s) of (jl.ip . which is real for real s and satisfies 

q(s) = Ai(s)(l + o(l)), as s — ^ -t-oo, (1.2) 

where Ai denotes the usual Airy function. 

The Hastings-McLeod solution q(s) (or functions related to it) do not appear in the local 
scaling limits of the correlation kernel at time t ^ t CI a. The Hastings-McLeod solution, however, 
does appear in the asymptotics of the recurrence coefficients of the multiple Hcrmite polynomials 
related to the non-intersecting Brownian motions, see [19 . 

The results in [T!5] were obtained from the Deift-Zhou steepest descent analysis of the 4x4 
matrix-valued Riemann-Hilbert (RH) problem for multiple Hermite polynomials. During the 
analysis we had to construct a local parametrix at a point xq € K that lies strictly between 
the two intervals where the two groups of Brownian paths accumulate. The point xq does not 
have a physical meaning. However, the local parametrix affects the recurrence coefficients of the 
multiple Hermite polynomials, as mentioned before. 

The aim of this paper is to perform a similar steepest descent analysis for the critical time 
t = icrit- This multicritical situation, after appropriate scaling of the parameters will be locally 
described by the solution of a model RH problem of size 4x4. The RH problem can be considered 
as a combination of two RH problems for the Airy function with an additional non-trivial coupling 
in the jump matrices. 

Using this new 4x4 model RH problem, we obtain an expression for the limit of the correlation 
kernel at the critical time. We show that the kernel has an integrable form determined by the 
solution to the RH problem, see Theorem 12.71 

We find it remarkable that this new 4x4 RH problem is again related to the Hastings-McLeod 
solution of the Painleve II equation. More precisely, we prove that the Hastings-McLeod solution 
shows up in the residue matrix in the asymptotic series at infinity of the 4x4 model RH problem, 
see Theorem 12.41 This is very similar to the situation for the classical 2x2 RH problem for 
Painleve II due to Flaschka- Newell [HI [53]. This suggests that our 4x4 problem may be 
expressible in terms of this smaller 2x2 problem; however we could not find such an expression. 
So our RH problem might lead to a genuinely new Lax pair for the Hastings-McLeod solution to 
Painleve II. See also the paper [30] for Lax pairs with matrices of size 3x3. 

As a consequence, we are able to show that the results in [19] on the recurrence coefficients 
remain valid at the critical time t = i cr ;t, i.e., the asymptotic behavior of the recurrence coeffi- 
cients of the multiple Hermite polynomials is still governed by the Hastings-McLeod solution to 
Painleve II with exactly the same formulas as in |19j . We find this a surprising fact. 

Very recently, a model of non-intersecting random walks was studied by Adler, Ferrari and 
Van Moerbeke [4] in a situation that is very similar to ours. There are two groups or random 
walks in [4], that in the scaling limit fill out two domains that are tangent in one point. It is 
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shown that there is a limiting correlation kernel at the tacnode, and two expressions for it are 
given in terms of multiple integrals involving Airy functions and the Airy kernel resolvent. It 
seems very likely that this limiting kernel should be equivalent to the one that we obtain in 
Theorem 12 . 71 below for the symmetric case (i.e., p* = p\ = 1/2 in Theorem 12. 7p . but we have not 
been able to make this identification. 

It would indeed be interesting to see how the Painleve II equation arises in the framework of 
[3], and conversely, to see how our formula can be reduced to integrals with Airy functions and 
the Airy kernel resolvent. It would also be interesting to have a process version with an extended 
tacnode kernel. 



2 Statement of results 

We now give a precise statement of our results. In Sections 12.1112.31 we describe our situation 
and the connection with a 4 x 4 matrix valued RH problem. In Section 12.41 we formulate the 
the new 4x4 RH problem and its properties, in particular the connection with the Painleve II 
equation. Finally, in Sections I2.5H2.6I we state the main results about the limiting behavior of 
the correlation kernels and recurrence coefficients. 



2.1 Correlation kernel and the Riemann-Hilbert problem 

We consider n one-dimensional non-intersecting Brownian motions with two starting points a\ > 
a 2 at time t — and two ending points b\ > b 2 at time t — 1. We assume that ni of the particles 
move from the topmost starting point a\ to the topmost ending point b\, and n-i particles move 
from the bottommost starting point a 2 to the bottommost ending point b 2 , with n\ + n 2 = n. 
The transition probability density of the Brownian motions is 

Pfr(t,x,y)= JL exp(--± T (x-y) 2 ) , a% = l/N, (2.1) 

with an overall variance a 2 N — 1/N, that decreases as n increases such that 

T = na 2 N = 1 > (2.2) 

remains hxed. We interpret T as a temperature variable. 

Using the above setting, it is known that the positions xi, . . . , x n of the Brownian paths at 
time t £ (0, 1) have a joint probability density [31] 

p(x u . ..,&») = -J- det (fi(xj))" j=1 det {g%{xj))^ j=1 (2.3) 



with functions 



gi-l 

fi(x) = gxil P N i.t,a 1 ,x), i=l,...,m, (2.4) 
gi-l 

fm+i(x) = gxi _ 1 P N (t,a 2 ,x), i = l,...,n 2 , (2.5) 
gi-l 

9i{x) = dxi _ x PnO- ~ t, x, h), i = l,...,m, (2.6) 

gi-l 

ffm-NO) = q i-x P N(l ~ t,x,b 2 ), i = l,...,n 2 , (2.7) 
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and with Z n a normalization constant. Note that (|2.3[) is a biorthogonal ensemble [10]. In 
particular, it is a determinantal point process with correlation kernel 



K n (x,y)= ]T f i {x){A- 1 ). .g j {y), (2. 



where (A -1 )^. denotes the (i,j)th entry of the inverse of the matrix 

A = (/ ^ X ^ 9j ^ dx ) 

The double sum formula (|2.8p for the kernel is not very tractable for asymptotic analysis. 
However there is a convenient representation of K n in terms of the solution of a Riemann-Hilbert 
problem. It is of size 4x4, since 4 is the total number of starting and ending positions. 

Define the weight functions 

w hk {x) =exp(-^(x 2 -2 afc2; )) , fc=l,2, (2.9) 



wa,<(ar) = exp I - _ (x^ - 26^)1 , J = 1,2, (2.10) 

and consider the following RH problem which was introduced in |14) as a generalization of the 
RH problem for orthogonal polynomials [24], see also [45] . 

RH problem 2.1. We look for a 4 x 4 matrix-valued function Y : C \ K — > C 4x4 satisfying 

(1) i^(z) «s analytic for z € C \ K. 

(2) y /ias limiting values Y± on K, where Y+ (Y-) denotes the limiting value from the upper 
(lower) half-plane, and 

Y + (x)=Y_(x)^ xeR, (2.11) 

where 1% denotes the 2x2 identity matrix, and W(x) is the rank-one matrix (outer product 
of two vectors ) 

W(x) = (^^) Ki(tf) W2 , 2 (x)). (2.12) 

(3) .As z — > oo, we /law £/ia£ 

y(«)= ^/+^ + 0^^diag(z Ill ,z n2 ,z-" 1 ,z- n2 ). (2.13) 

The RH problem has a unique solution that can be described in terms of certain multiple 
orthogonal polynomials (actually multiple Hermite polynomials of mixed type), see [TH [19] for 
details. According to [T3], the correlation kernel K n from (|2.8p is equal to 

( toi,i(a:)^ 
101,2 (a:) 


V o / 

The kernel K n in (|2.14l) not only depends on n, but also on n\, n-i, a\, 02, b±, &2, t and T. In 
what follows, we will take the variables T and t fixed while the other variables will be varying 
with n, see Section! 



K n {x,y) = 2ni( l_ y) (° «»Av) w 2 Ay))Y+ 1 (y)Y + {x) 



(2.14) 
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2.2 Separation of the starting and ending points 

We now discuss in more detail the three situations in Figure [1] Given n\ and n 2 , we denote the 
corresponding fractions of particles by 

Pi-= — , P2 ■= — =l-pi, (2.15) 
n n 

which are varying with n, and we assume that 

Pl =pl+0(l/n), P2 =p* 2 +0{l/n), n->oc, (2.16) 

for certain limiting values Pi,p 2 € (0, 1). Of course, p 2 = 1 — p\. 

For given T and pl,p 2 , the three cases of large, small and critical separation of the starting 
and ending points are distinguished as follows; see [19]. There is large separation in case the 
inequality 

(01 - 02)(6i - 62) >T(^* + ^|) 2 (2.17) 

holds. Then the Brownian motion paths remain in two separate groups, and the limiting hull in 
the to-plane consists of two ellipses. For any t € (0, 1), the limiting distribution of the positions 
of the paths at time t is supported on the two disjoint intervals [aJ,/3*] and [aj,/^], where the 
endpoints satisfy 

a* 2 < ft < a\ < PI (2.18) 

and are given explicitly by 



a* = (l- t)a 3 + tbj - 2Jp*Tt(l - t), (2.19) 



/3* = (1 - t)a 3 + tbj + 2^Jp*Tt(l - t), (2.20) 
with the limiting density on these intervals given by the semicircle laws 

2l <Tt{l-t) ^l { ^ {t) ~ x){x ~ a *i {t)) ' (2.21) 

for j = 1,2. See again Figure [TJa). 

We are in the case of small separation if the inequality (|2.17l) is reversed, so that 

(ai-a 2 )(h-b 2 ) <t(vpI + VpI) 2 ■ ( 2 - 22 ) 

Then the limiting hull in the te-plane is an algebraic curve and the limiting distribution of the 
paths at any given time t 6 (0, 1) is not given by semicircle laws anymore. This was shown in 
[IB] for the case where p\ = p 2 = 1/2. See Figure Hfb). 

The case of critical separation corresponds to the equality sign in (|2.17[) . that is, 



(01-02X61-62) =t(vp! + VpA) 2 ■ (2- 23 ) 



Then the limiting distribution is still given by two semicircle densities (|2.21[) . which are now 
however touching at the critical time t CT i t £ (0, 1) defined by 

^ = 1 — a \Td h y ( 2 - 24 ) 

(ai - a 2 ) + (01 - b 2 ) 
in the sense that j3 2 = a \ in (|2.18l) when t — t CT i t . See Figure QJc). 
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2.3 Double scaling limit 

In order to study the case of critical separation we choose without loss of generality 

T = l, (2.25) 

and take a\ > a 2 , > b 2 so that 



cf. (ED , We also take 



and 



(a* 1 -a* 2 )(b* 1 -b* 2 ) = [VpI + V^)\ (2-26) 



t - tait ~ (al-a*) + (lt-b* 2 ) (2 ' 27) 

£crit = (1 — icritj ^ V ^crit ~ (2.28) 



as the critical time and place of tangency. 

We can assume that T and t remain fixed as in (12.25)) and (|2.2Z[) without loss of generality, 
since we cannot create different limiting behavior by varying T and t as well, see Remarks 12. Ill 
and l2~T2l in Section [231 

We let the starting points oi, a 2 and ending points b±, b 2 be varying with n in such a way 
that for certain real constants L\, L 2l £3, L4, 

01 = a{ + L in - 2/3 , b x = b\ + L 3 n" 2/3 , (2.29) 

a 2 = a* 2 + L 2 n- 2/3 , b 2 = b 2 + L 4 n" 2/3 , (2.30) 

as n — > 00. It turns out that the constants only appear in our results through the combinations 

L 5 := (ft? - b* 2 )L x + (aj - a* 2 )L 3 , (2.31) 
L 6 := (bl - b* 2 )L 2 + (oj - a* 2 )L 4 . (2.32) 

Near the critical value x C rit, the asymptotics of K n for n — > 00 will be described by a family 
of limiting kernels 

K tacnode {u,v;r 1 ,r2,a 1 ,a a ), u,veR, 

depending on four variables 7*1, r 2 > and si,s 2 € K that depend on the values p*, -^5 an( i 
Because of a dilation and translation symmetry 

c 2 K tacnode (c 2 u, c 2 v- n, r 2 , si, s 2 ) = X tac " ode ( M , w; c 3 ri , c 3 r 2 , csi, cs 2 ), c > 0, 

K tacnode( u + 2Ci v + 2c; n,r 2 , s 2 ) = K tacnode (u, v; r u r 2 , si - cn,s 2 + era), c € K, 

the family essentially depends on two parameters only, and we could for example choose r± = 1 , 
si = s 2 . However, in order to preserve the symmetry in the formulas, we prefer to use four 
parameters. 

The limiting kernels are given in terms of a RH problem that we discuss next. 
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Figure 2: The figure shows the jump contours in the complex £-plane and the corresponding 
jump matrices Jk, k = 0, . . . , 9, in the RH problem for M = M{Q. 



2.4 A new 4x4 Riemann-Hilbert problem 

The 4x4 matrix- valued RH problem has jumps on a contour in the complex plane consisting of 
10 rays emanating from the origin. The rays are determined by two numbers ipi, tp2 such that 

< if! < ip 2 < tt/3. (2.33) 

The value 7r/3 in (|2.33p is not the optimal one, but it will be sufficient for our purposes. We 
define the half-lines Tk, k = 0, . . . , 9, by 

r =M+, T 1 =e 1Vl R + , T 2 = e lV2 R + , 

■t \ -l ■/ . (2.34) 

and 

r 5+fc = -r fc , fc = o,...,4. (2.35) 

All rays are oriented towards infinity, as shown in Figure [2] 

RH problem 2.2. We look for a 4 x 4 matrix-valued function Af(£) (which also depends para- 
metrically on the parameters ri,r% > and Si,sa G C) satisfying 

(1) M(C) ts analytic for C G C \ (uLo r fc) ■ 

(2) For ( £ Tt, the limiting values 

M+(C)= lim M(z), M_(C)= lim M(z) 

2:— ^C,^ on +-side ofT^ z—>i^,z on — side ofT^ 
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exist, where the +-side and — -side ofTk are the sides which lie on the left and right ofTk, 
respectively, when traversing Tk according to its orientation. These limiting values satisfy 
the jump relation 



M+(C)=M_(c)J fe) Cer fe , k = o, 

where the jump matrices Ju are shown in Figure 
(3) As ( — > oo, we have that 



Q 



(2.36) 



V2 



/ 1 -i 

1 i 

-i 1 

\0 i 1/ 



jdia g ((-c)- 1/4 ,r 1/4 ,(-c) 1/4 ,c 1/4 ) 

dia^e"* 1 ^, e^O^O e 9 ^), (2.37) 



where the coefficient matrices Mi, Ma,... depend on the parameters r\, r 2 , s±, S2, but not 
on C, and where we define 



and 



(4) M(C) is bounded near £ = 0. 



l(C) = ^ 1 (-C) 3 / 2 + 2 Si (-c) 1 / 2 



9 2 (C) = §r 2 C 3/2 + 2 S2 C 1/2 - 



(2.38) 
(2.39) 



In (I2.37p ^ (|2.39p . we use the principal branches of the fractional powers, so that for example 
(— C) 1 / 4 is defined and analytic for £ G C \ [0, oo) with real and positive values on R~. We write 

M(C;ri,r 2 ,si,S2)) 

in case we want to emphasize the dependence of M on the parameters. 

It follows from standard arguments (e.g. [TB]) that the solution to the RH problem 12.21 is 
unique if it exists. The existence issue is our first main theorem. 

Theorem 2.3. (Existence:) Assume that ri,ra > and Si,sa € M. Then the RH vroblem 12.21 
for M(C) is uniquely solvable. 

Theorem 1 2 . 31 will be proved in Section [4] by using the technique of a vanishing lemma (171 1251 



Our next result is about the residue matrix M\ in (|2.37p . We show that its top right 2x2 block 
is related to the Hastings-McLeod solution q(s) of the Painleve II equation and the associated 
Hamiltonian u(s). 

Theorem 2.4. (Mi vs. the Painleve II equation:) Let the parameters ri,ra > and S\,S2 € K 
in (|2.37[) - (|2.39[) be fixed. Then the matrix M\ in (|2.37[) can be written as 



Mi 



a b ic id\ 

—b —a id ic 

ie if —a b 

\if ie —b a/ 



(2.40) 
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for certain real valued numbers a, a, b, b, c, c, d, e, e, / that depend on n, T2, Si, #2- ^ n addition, we 
have 



(rir 2 )V6 

2/3 



s 2 

u(tr) + -i, (2.42) 



rj /3 (rf + r|)V3 n 

2/3 



u{a) + (2.43) 



r2 /3 (r 2 +r|)i/3 r 2 

where q(s) is the Hastings- McLeod solution of the Painleve II equation (ll.lj) - (|1.2l) . «(s) is the 
Hamiltonian 

u( S ):=(q'(s)) 2 -sq 2 ( S )-q 4 (s) (2.44) 

and 

2(ri so + rosi ) 

(2.45) 



' (rirajVa^ + ^jVs- 

Theorem 12.41 shows that the top right 2x2 block in (|2.40p behaves in a similar way as the 
residue matrix in the classical 2x2 matrix-valued RH problem for Painleve II due to Flaschka 
and Newell [22l [23] . We do not think that our 4x4 RH problem can be reduced to the 2x2 
problem. 

The proof of Theorem 12.41 is given in Section [5] 

Remark 2.5. In Section [5] we also find identities relating the entries in the top left block of Mi 
to the entries in the top right block in (|2.40l) . namely 

dd v 

(2a + c 2 )r l = r 2 d 2 + Sl , (2a + c 2 )r 2 = nd 2 + s 2 , — = 2cd - 2b = —(2cd - 2b) 

OSi 7*1 

see (|5.9p . (|5.10[) , and (|5.22p . Hence also a,a,b and b can be directly expressed in terms of q(s). 
We also obtain the Lax pair equations 



DM _ uM dM 
dQ dsj 



UM, — = V j M, j = 1,2 (2.46) 



with matrices U and Vj that are explicitly given in terms of Mi, see Propositions 15.31 and 15.51 In 
fact the matrices only depend on the entries a, a, b, b, c,c and d of Mi, and so U and Vj depend 
only on the Hastings-McLeod solution q(s) of Painleve II 

, 

asymptotic behavior (|2.37[) in that sector. 



In terms of U one may therefore view M as a solution of ^M- = UM in each sector with the 



2.5 Critical limit of correlation kernel 

The limiting kernels are defined in terms of the solution of the model RH problcm l2.2l as follows. 
Let M(C; ri, r2, si, S2) be the solution of the RH problem 12.21 for fixed ri,r2 > and si,s 2 € M. 
Thus, M is analytic in each of the sectors determined by the contours T^, and the restriction 
to one such sector has an analytic continuation to the entire complex £ plane. In other words, 
the entries of M are entire functions. Consider the sector around the positive imaginary axis, 
bounded by T2 and we denote the analytic continuation of the restriction of M to this sector 
by M. 
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Definition 2.6. For u,«el, the kernel K tacn (u, v; ri, r2, si, S2) is defined by 



K 



tacnode 



(u,v;n,r 2 ,si,s 2 ) 



1 



2ni(u — v) 



(0 1 1) M 1 (u;ri,r 2 ,Si,s 2 )M(v;ri,r 2 ,Si,s 2 ) 



1 





(2.47) 



We can rewrite the kernel in terms of the limiting values M + (u) and M + (v) of M on the real 
axis, by using the jump relations in the RH problem for M . For example, for u, v > 0, we have 



K 



tacnode 



(u,v;r 1 ,r 2 ,s 1 ,s 2 ) 



2-ki{u — v) 



(-10 1 0) Af+ 1 (u)Af+(w) 




1 



(2.48) 



with a different expression in case u and/or v are negative. Here, we have dropped the dependence 
of M on the parameters 7*1, r 2 , Si, s 2 . 

From (|2.47l) , it is also easily seen that the kernel has the integrable form 

h(u)gi{v) + f2(u)g2(v) + f3{u)g 3 (v) + fi(u)g<i(v) 



for certain entire functions fj and gj, j = 1, . . . ,4, with X)j=i fj( u )9j( u ) — 0- The functions of 
course depend on ri,r 2 , s%, s 2 . 

The following is the main theorem of this paper. 

Theorem 2.7. (Correlation kernel at the tacnode:) Consider n non-intersecting Brownian mo- 
tions on [0, 1] with transition probability density (|2.1[) with two given starting points a\ > 02 and 
two given endpoints b\ > b 2 . Suppose n\ paths start in a\ and end in b\, and n 2 = n — n\ paths 
start in a 2 and end in b 2 - Assume that T = 1 and that (|2.16p . (I2.29[) - (|2.30p hold with values 
pl,p 2 , a\, 02, &!, &2 such that (I2.26[) holds. 

Let t cr it and x cr u be the critical time and place of tangency as given in (|2.27l) and (|2.28D . 
respectively. Then the correlation kernels K n for the positions of the paths at the critical time 
t C rit satisfy 



lim 

n— >oc 



exp (c2n 1 ^ 3 (it — «)) 



7 2 /3 



/ U V \ 

I ^ crit T" TTT , % crit T" q~7T I 

V cn 1 J cn^i^J 



= K 



tacnode 



(u,v;ri,r 2 ,si,s 2 ), (2.49) 



with K tacnode given by (12T4T1) and 

I *\l/4 

n = (Pi) , 



a/4 



si = 



2 (\ 






/p! + 




and 


C2 m 


(|2.49|) c 




1 





r 2 = (^) 1/4 , 

(P2) 1/4 



S2 



2(\M + v/pI) 



\/t critic tcrit) 

c 2 := c 



/1 \ a i &i 

(1 <>crit) tcrit~^ 



0*2 , , ^2 



(l tcrit) 2 ~^ t cr it ^ 



(2.50) 
(2.51) 

(2.52) 
(2.53) 
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Theorem 12 . 71 will be proved in Section I8TT1 

Remark 2.8. The extra factor exp(c2?i 1 ^ 3 (u — v)) in (|2.49[) is irrelevant as it does not change any 
of the determinantal correlation functions det (K n (xi, Xj)) of the determinanal point process. 

Remark 2.9. In view of ([235]) . (|23T1) and ([2T3"2"]l . we may rewrite (j2~5T]) as 





- a* 2 


K 


-K 



31 - y ( y \^ 2 Ll + V t^f 7 ' 52 " "? ( V ^Mt i2 + 1 / ^^ L4 

Remark 2.10. With the values (|2.50l) ~ (|2.5ip . the parameter cr of (|2.45|) that is the argument of 
the Painleve II transcendent is related to p\, p 2 , L§ and Lq by 

(P ^ )1/6 (iB-ie)- (2-54) 



(VpJ + v^l) 4/3 

Remark 2.11. Instead of varying the endpoints, we could have used the temperature T as a 
scaling parameter with the same scale 0(n~ 2 / 3 ). If we keep a%, a 2 ,b\,b 2 fixed at the values 
, 03 , 6J , 62 an( i vary T with n such that 

T = 1 - L 7 n- 2/3 +o(n~ 2/3 ) as n -> 00, 

then it can be checked that the same local behavior at the tacnode can be created with a change 
in the endpoints while keeping the temperature fixed. Indeed, to that end, we vary the endpoints 
as in P^g]) - ([2T50|1 with 

Li = ^a\L 7 , L 2 = ^a 2 L 7 , L 3 = l a = \ h 2 L i- ( 2 -55) 

Thus, by ([2~3"Tj) and (j2~32l) . it follows that 

L 5 = ^((bl-b* 2 )a* 1 + (at-a* 2 )bt)L 7 , 

L 6 = \{{bl-b* 2 )a* 2 + {al-a* 2 )bl)L 7 , 

and 
so that 

^ = (piP2) 1/6 (v^i"+v^) 2/3 ^, 

with the aid of (12.541) . This expression for a is compatible with the one in [T!5] . 

Remark 2.12. We could similarly vary the time t around the critical time t cr it while keeping 
other parameters fixed, say, 

t = t crit + L s n- 2/3 + o(n- 2/3 ). 

This effect can be modeled by varying the endpoints as in (|2.29l) - (|2.30p with t = t cr j t and T = 1 
fixed, with now 

I a i j r b l r 

Ll ~2t crit (i-t crit ) 8 ' 3 ~ ^Fu) 8 ' 

a* 2 b* 2 (2 - 56) 

-^2 = — ~T TZ T-^S, -^4 = 7^ T\ 7-^8- 

2^crit(l ^crit ) 2£ cr jt(l ^crit ) 

Hence L^—Lq = and so cr = by (I2.54[) . Thus the argument (|2.45|) of the Painleve transcendent 
does not change if we only vary t. 



13 



2.6 Additional results 



2.6.1 Semicircle densities 

There are two other results that come out of our asymptotic analysis. The first one is about the 
limiting density of the non-intersecting Brownian paths. 

Theorem 2.13. (Touching semicircle densities:) Consider the double scaling limit as described 
in Section \2.Sl Then as n —¥ oo, the Brownian particles at the critical time t in (|2.27[) are 
asymptotically supported on the two touching intervals [o^,/?^] and [a^,/^] fl2.19[l - (|2.20p . with 
02 = a* , and with limiting densities given by the semicircles (|2.21[) . 

Theorem l2.13l can be proved from the steepest descent analysis in Section [7] in quite the same 
way as in [19] . We will not go into the details. 



2.6.2 Critical limit of recurrence coefficients 

The second result concerns the recurrence coefficients of the multiple Hermite polynomials. In 
the critical limit, their behavior is governed by the Hastings-McLeod solution of the Painleve II 
equation, in exactly the same way as in [19] . For the 4x4 matrix Y\ in (|2.13|) , the combinations 
(Yi)ij(Yi)j j i, i, j = 1, ... ,4 with i ^ j are called the 'off-diagonal recurrence coefficients' in [19] . 
where we use the notation (M),-^ to denote the (i,j)th entry of any given matrix M. It is also 
shown in [TpJ that the quantities (Yi)i,j(Yi)j,i w ^ tn ? - 3 are determined by (^1)1,2(^1)2,1 and 
(^1)1,4(^1)4,1 alone. We then have the following result. 

Theorem 2.14. (Asymptotics of off- diagonal recurrence coefficients:) Consider the double scal- 
ing limit as described in Section \2.'6\ Then we have 

(^i) 1,2(^1)2,1 = -KHl rit {b\ - blfq\a)n-^ + O^ 1 ), (2.57) 
(^1)1,4(^)4,1 = K 2 t mt (l - t crit )(al - a* 2 )(bl - b* 2 )q 2 (a)n- 2 / 3 + 0(0, (2.58) 

as n — >• 00, where q(s) is the Hastings-McLeod solution to the Painleve II equation, a is given by 
([2754]) and 

K = . (2 . 59) 



Theorem 12.141 will be proved in Section [ 

Remark 2.15. For ease of comparison, we have stated the above theorem in exactly the same 
way as in [19]. We note however that the expressions (I2.57P and (|2.58l) can be simplified in the 
present case, since we are now looking exactly at the critical time t — t cr - lt . This means that 
the variable t in the above formulas can be substituted by (I2.27[) . Further simplification can be 
obtained by substituting (|2.26[) . 

Remark 2.16. In [T!5] there is also a result on the Painleve II asymptotics of the so-called 'diagonal 
recurrence coefficients'. One can show that exactly the same result holds in the present case. 
This will be briefly commented at the end of Section [ 



2.7 About the proofs 

The rest of this paper contains the proofs of the theorems and is organized in two parts: Part I 
(Sections and Part II (Sections [6] -EJ}. Part I deals with the RH problem [221 for M{Q. It 
consists of three sections: In Section [3] we perform an asymptotic analysis of this RH problem 
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for s\ —¥ +00; in Section 0] we establish the existence result in Theorem l2.31 and in Section [5] we 
prove the connection with the Painleve II equation in Theorem 12.41 

Part II deals with the critical asymptotics of the non-intersecting Brownian motions at the 
tacnode. The proofs of Theorems 12 . 71 and 12 . 141 will be given in Section [51 They are based on the 
Deift-Zhou steepest descent analysis of the RH problem l2.il which will be discussed in Section[7] 
In the proofs an important role is played by two modified equilibrium problems that give rise to 
two functions £1, £2 and their anti-derivatives Ai, A2. The modified equilibrium problems will be 
discussed in Sectional 

The two parts are largely independent. Both parts contain the steepest descent analysis of a 
4x4 matrix valued RH problem, that we give in some detail, although it makes the paper rather 
lengthy. Some of our notation will have different meanings in the two parts. For example, V\ 
and V 2 are used to denote certain 4x4 matrices in part I, see (|5.14l) and (|5.15p . while in part II 
they denote two external helds in an equilibrium problem, see (|6.ip . We trust that this will not 
lead to any confusion. 

Part I 

Analysis of the Riemann-Hilbert 
problem for M(£) 

3 Asymptotic analysis of M(£) for s\ — > +00 

In this section we analyze the model RH problem l2.2l for M(() as s\ — > +00. Our goal is twofold: 
we want to prove the solvability of the RH problem for si sufficiently large, and establish the 
large si asymptotics for the quantities c, c and d in (|2.40l) . The goal of this section is to prove 
the following proposition. 

Proposition 3.1. Let r%,r2 > and s 2 <E R be fixed. Then for large enough s± £ R, the RH 
vroblem \2.2\ for M(C; r1.r9.s1. sn) is uniquely solvable. 
Moreover, we have 

d := -i (Mi(ri, r 2 , si, s 2 ))i j4 

^(^) V4 -(-H (i+o(a ' 3/S)) - (3 ' 1) 

s 2 

c := -i (Mi(n,r 2 , ai, s 2 )) 1 3 = -± + O (s^ 1 ) , (3.2) 

n 
s 2 

c:= -z(M 1 (r 1 ,r 2 ,s 1 ,s 2 )) 24 = -J + (9(sf 1 ), (3.3) 

os si +00 , where a is given by (|2.45p . 

We will use Proposition 13 . 1 1 in the proofs of Theorems 12.31 and 12.41 

In the proof of Proposition 13.11 we will hrst show that we can restrict ourselves to the case 

n = 1, s 2 = 0. (3.4) 

This is due to certain dilation and translation symmetries for M. The further analysis will then 
be based on the Deift-Zhou steepest descent method, using ideas of [THIHZ]- For the problem at 
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hand, the analysis consists of a series of transformations M i->- A n- B M> C M> D, so that the 
matrix- valued function D uniformly tends to the identity matrix as si — > +00. 

Remark 3.2. Although we will perform the transformations for s\ > and sufficiently large in 
what follows, it is also possible to apply an asymptotic analysis as s± — ¥ —00. But in that case 
the 'g-functions' will be more complicated: They must be constructed from a 4-sheeted Riemann 
surface of genus zero, which has a somewhat similar flavor as the one in [IS]. We will not go into 
the details. 



3.1 Reduction to the case r\ = 1 and s 2 = 



The four parameters rx,r 2 , si, s 2 in the RH problem 12.21 for M are not independent. Indeed, we 
have the following dilation symmetry 



/ 7 l/2 






















T -l/2 


V 









\ 




-1/2 J 



M(7 2 C;ri,r 2 ,si,s 2 ) 



^(C;7 3 n,7 3 ^,7Si,7S2), 7>0, (3.5) 



and the translational symmetry 



1 


i8(2s\ — r\S) 




0\ 



1 



iS(2s 2 + r 2 S) I) 



M(C + 2(5;r 1 ,r 2 ,s 1 ,s 2 ) 

w M(C; ri,r 2 , s\ - r%6, s 2 + r 2 S), 



(3.6) 



where the symbol ~ in (|3.6[) is used in the sense of 'equality up to contour deformation'. The 
equality (13.51) is immediate, due to the fact that both sides satisfy the same RH problem. To 
show p. 6)1 . we note that the left-hand side of (13.6)1 satisfies a RH problem with the same jumps 
as M, but on a shifted contour. By an easy transformation of the RH problem, we may shift the 
contour back to the original contour. After this transformation, the equality Q3.6P holds, since 
both sides of (|3.6t also have the same asymptotic behavior as C ~> 00. 

In view of Q3.5P and Q3.6)) . we may restrict ourselves to r\ — 1 and s 2 = in order to prove 
the solvability statement in Proposition 13. II 

For the residue matrix Mi in ()2.37|i . we find from (13.51) that 

Mi(7 3 ri,7 3 r 2 ,7Si,7s 2 ) = 

-i diag (7 1 / 2 , 7 1/2 , 7~ 1/2 , 1' 112 ) M 1 (r 1 ,r 2 ,s 1 ,s 2 ) diag (7-1/2, 7 -i/a, 7 V2, 7 V2) , (3 . 7) 



16 



and from (13. 6[) . after straightforward calculations, 



Mi(ri,r 2 , si - ri<5, s 2 + ?" 2 <5) = Mi ,5 

/ 1 0^ 

10 

i5(2s 1 -r 1 S) 1 



V 







Mi(ri,r 2 ,si,s 2 ) 



iS(2s 2 + r 2 5) 1/ 



1 



-iS(2 Sl - r x S) 




0\ 



1 



V 



-iS{2s 2 + r 2 S) 1/ 



where 



(Mu) M = -(M M ) 3 , 3 = - 2s 2 6 2 + 2r lSl 6 3 - ^r\5\ 

(M M ) 2 , 2 - -(M M ) 4 ,4 = + 2s^ 2 + 2r 2 s 2 5 3 + X -r\b\ 

(-^1^)1,3 = i(—2si5 + riS 2 ), 
(M 1 j) 2 , i ^i{2s 2 S + r 2 S 2 ), 

(Mi, s ) 3 ,i = i f- Sl 5 2 + i(2n - 8s 3 )5 3 + 4 s 2 r 1( 5 4 _ 2sir 2 5 s 



(3.8) 



(M M ) 4 , 2 = i \yS2S 2 + g(-2r- 2 + 8s^)(5 3 + 4r 2 s 2 5 4 + 2r 2 s 2 S 5 + -r^ 6 

and the other entries of M± t s are zero. 

In view of the structure of Mi in (12.40[) . a combination of (|3.7I) and (|3.8p yields 

1 

d(7 n,7 r 2 ,7(si - n<5),7( s 2 + r 2 d)) = -d(n,r 2 ,si,s 2 ), 

7 

c(7 3 ri,7 3 r 2 ,7(si -n J), 7(52 + r 2 <5)) = - (c(ri, r 2 , si, s 2 ) - 2si<5 + riS 2 ) , 

c(7 3 ri,7 3 r 2 ,7(si -r 1 S),-f(s 2 + r 2 5)) = - (c(ri, r 2 , si, s 2 ) + 2s 2 <5 + r 2 <5 2 ) , 

for 7 > and <5 6 ML These relations are consistent with the formulas for d, c and c given in 
p.l[) - (|3.3p as well as those in (|2.41[) - ([2.45p . Thus, we may indeed restrict to the case r\ = 1 and 
s 2 = in the proof of Proposition l3.il 



3.2 First transformation: M ^ A 



We consider the RH problem 12.21 for M with parameters r\ = 1 and s 2 = 0. The first transfor- 
mation is a rescaling of the RH problem. 



Definition 3.3. We define 



A(0 = dia g ( s ; /4 , s y\ S " 1/4 , s - 1/4 )M( Sl c), C e c \ ( |J r fc j 

\fe=0 / 



(3.9) 



Then A satisfies the following RH problem similar to that for M. 
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RH problem 3.4. 

(1) is analytic for £ 6 C \ ^lJfc=o r^J , where Tk, k = 1, 9 is shown in Figure^ 

(2) A /ias i/ie same jump matrix Jk on as M . 

(3) ^4s £ — > oo, we Ziave 





r 1 





—i 


°\ 


1 





1 





j 




— i 





1 









i 





V 



with 
and 



l))diag((-o- 1 / 4 ,r 1/4 ,(-o 1/4 ,c 1/4 ) 

diag(e- (0 , e - Ae " a (C) , e** 1 (C) , e Ae ~ 2 (C) ) , (3.10) 

A = s 3 / 2 , (3.11) 

*i(0 = §(-C) 3/a + 2(-C) 1/2 , 02(0 - ^r 2 C 3 / 2 . (3.12) 

The residue matrix Aj in (|3.10l) is given by 

Ax = — diag(s} /4 , s} /4 , s^ 1/4 , s7 1/4 )Mi diag(s7 1/4 , s^ 1 ^, s} /4 , s} /4 ). (3.13) 
si 

In view of (|3.13[) , we find the following expressions for the numbers c, c and d in (|2.40p : 

c = -iV^i(^i)i,3) c = -iy/sl(Ai) 2 ,i, d = -iv^i(^i)i,4- (3-14) 

3.3 Second transformation: 

In the second transformation we apply contour deformations; see also [27 . The six rays Tk, 
k = 1, 2, 3, 7, 8, 9, emanating from the origin are replaced by their parallel lines emanating from 
some special points on the real line. More precisely, we replace Ti and Tg by their parallel 
rays Ti and Tg emanating from the point 2, replace T 2 and T% by their parallel rays T 2 and Tg 
emanating from the point 2 — e, and T3 and I\ by their parallel rays I^ and I\ emanating from 
e, where e is a small positive number. See Figure [3] 

Definition 3.5. Denoting by Eij the 4x4 elementary matrix with entry 1 at the (i, j)th position 
and all other entries equal to zero, we then successively define 

A(()(I + E 3A ), for C between i\ and fi, 
A(£)(I ~ Es A ), for C between T 9 and f 9 , (3.15) 
A((), elsewhere, 

B^(()(I -E2.1+ E 3A ), for C between T 2 and f 2 

and for ( between T 8 and T$, (3.16) 
B^'(Q, elsewhere, 




and 



(C) {I + £1,2 - £4,3) > for C between T 3 and f 3 
B(C) = { and for ( between T 7 and f 7, (3.17) 

i?( 2 )(C), elsewhere. 
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Figure 3: Contour E^ for the RH problem 13.61 for B. Note the reversion of the orientation of 
some of the rays. 



It is easily seen that B is analytic in C \ E# , where E^ is the contour shown in Figure |3] Note 
that we reverse the orientation on some of the rays, in particular the real line is oriented from 
left to right; compare Figure [3] with Figure[2] Moreover, B satisfies the following RH problem. 

RH problem 3.6. 

(1) -B(C) is analytic for ( e C \ E^. 

(2) B has the following jumps on E^: 



where Jb is defined by 



s + (C) = B_(C)Jfl(0 J 



MO = 

MO = 
MO = 
MO = 

MO = 

MO = 
MO = 
MO = 

MO = 



o 



-l 



(I + E 3 ,i), 
(I + E lt3 ), 

(I - E 2A i 

fl 1 

1 1 

1 

\0 

(I + E h2 - 
{I + E 2A ), 






V 



for C € (2, +00), 



for c g r x u r 9 , 

for CG (2-e,2), 
£3,4), for C G f 2 U f 8 

1 



V 

£4,3), /or C G f 3 U f 7, 
/or CG(0,e), 

for c g r 4 u r 6 , 



/orCG(e,2-e), 







\o -1 



1 





/or C G (-oo,0). 
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(3) As C — > oo, we /icwe 



s(0 



(l))dia g ((-c)- 1/4 ,r 1/4 ,(-c) 1/4 ,c 1/4 ) 



v/2 



/ 1 -i 0\ 

1 i 

-i 1 

\0 i 1/ 



diag(e- A91 (0 , e - ^ (0 , e ^ (0 _ e A* 2 ( C ) ) 



(3.18) 



where X — s 3 / 2 and 6i(£), $2(C) * s given in (|3.12p . 

By (|3.14p and the transformation A M- B in (|3.15D - (|3.17p . it is readily seen that 
c= -iy/si(Bi)x$, c = -i</si(Bi) 2 ,4, d = -i\/si{Bi)i^. 

3.4 Third transformation: B i— >■ C 



(3.19) 



In this transformation we partially normalize the RH problem 13.61 for B at infinity. For this 
purpose, we introduce the following '^-functions': 



and 



Note that 



si(C) = 3(2-C) 3/2 , CeC\[2,oo), 



52(C) = MO = g^C 372 , C e C\ (-00,0]. 



.91(C) = 0"i(O + (~Cr 1/2 - g(-C)- 3/2 + 0(C 5/2 ) 



(3.20) 

(3.21) 
(3.22) 



as C — > 00, where 9\ is given in (|3.12[) . 
Definition 3.7. We define 

C(0 = (I + i\E 3A )B(()dizg(e X!!l( -(\e X92( > c) ,e- Xgi ^ ,e- Xs2 ^). 
Then C satisfies the following RH problem. 
RH problem 3.8. 

(1) C(C) is analytic for ( £ C \ £3. 

(2) C has the following jumps on S^' 

C + (C) = C-(C)Jc(C), 



(3.23) 
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where Jc is given by 

MO 



MO 
MO 
MO 

MO 

MO 
MO 
MO 

MO 



/ 1 0\ 

10 

-10 0' 

\ 1/ 

(/ + e 2A Sl (C)£ 3l)j 

(I + e- 2X ^E h3 ), 



for(& (2,oo), 



for Ce (2-e,2), 
(Z + e^-ssXO^ - E 2;1 )), for C e f 2 U f 8 , 
/l e - 2A 9i(0 e -A(gi+g 3 )(C)\ 



5 -2A 92 (C) 



/orCe(e,2-e), 



/ 



1 e _A ^ fi,1+S2 ^'») 
1 
\0 

(Z + e -A(»-»)«}(£ lj2 - £ 4i3 )), / or C G f 3 U f 7 , 

(Z + e-^K)^), /or C 6(0,6), 

( / + e 2A S2(C ) S42); / or (er 4 ur 6 , 

/l 0\ 

1 

10 

\0 -1 0/ 



/or C e (-oo,0). 



(3) As £ — > oo , we /tawe 
C(C) 



Z- 



C 



o 



dia g ((-c)- 1/4 ,r 1/4 ,(-c) 1/4 ,c 1/4 ) 



V2 



( I -z 0^ 

1 

-i 

\ i 



i 

1 
1/ 

(3.24) 



where 



d = {I + iXE 3A ) 



Bi(I-iXEs,i) 



2 





-iX 


°\ 














a 

1 3 





A 2 
2 





V o 








0/ 



(3.25) 



Proof. The jump condition for C in item (2) follows from straightforward calculations, where we 
have made use of the facts that (<7i,+ + ,9i.-)(C) = for £ e [2, oo) and (52,+ + <?2.-)(C) = for 

Ce(-oo,o]. 

To establish the large £ behavior of C shown in item (3), we first observe from (|3.22[) and 

flnznn that 



^xonSim = 1 ± A (-C)- 1/2 + y(-C)- 1 T ^(-C)- 3/2 + 0(C 2 ), 



as £ — > 00, and 



e ±A( S2 -e 2 )(C) = 1 
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Hence, by (I3.18[) . it is readily seen that 



Affi(C) P A 92 (C) P -A S1 (C) P -A 32 (C)^ 



B(C)diag(e 



= ' / + f + ° (?)) dia s((-o- 1/4 ,r 1/4 ,(-o 1/4 ,c 1/4 )^ 



X 1 + 



A A 2 



A 



< 2C 3(-C) 3 / 2 



A A 2 



as C — » oo. Note that 



dia g ((-cr 1/4 ,c 1/4 ,(-c) 1/4 ,c 1/4 ) 



V2 



x / + 



A A 2 



-C 2C 3(-C) 3 / 2 



Ei, 



/=£ 2C 3(-C) 3/2 



/ 1 -i 0\ 

1 i 

-i 1 

\0 i 1/ 

A A 2 



/ 1 -« 0\ 

1 % 

-i 1 

\0 i 1/ 



3/2) 



-C 2C 3(-C) 3 / 2 



E; 



3.3 



/ 1 0\ 

10 

-iA 1 

y o o o 1/ 







-f 

V o o o o/ 



1 iA 



xdia g ((-c)- 1/4 ,r 1/4 ,(-o 1/4 ,c 1/4 )^ 



/ 1 -i 0\ 

1 i 

-i 1 

\0 i 1/ 



This, together with ([3~2"3) and (gUp) , implies (HOIl) and (JHHJ). 
It follows from (j3~T<?)) . (ET2"5|) and A = s 3/2 that 



(3.26) 



□ 



-iv/si(Ci)i, 3 + s 2 , c = -i v / si(C'i)2,4, d = 



(3.27) 



3.5 Asymptotic behavior of jump matrices for C 



The jump matrices of the RH problem 13.81 for C all tend to the identity matrix exponentially 
fast as A — > +oo, except for the jumps on (— oo,0) and (2, oo). This is easily seen for the jumps 
on f i U T 4 U (0, 2) U T 6 U f 9 , due to the facts that 

Rp „^J <0 ' if f <arg(C-2)<f , 

Ke 51 lU \ > 0, if < arg(C - 2) < f or f < arg(C - 2) < 2tt, l ^ 8j 



and 



i>2: ;;K?;r' <arsc<_1 ' <3 - 29) 



For the jump matrices on ^urs and TaU^, it is necessary to analyze the function (gi— <72)(C)- 
Recall that T 2 , r 8 (T3, T 7 ) intersect the real line at the point 2 — e (e, respectively) with e > 
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a small number. It is easily seen that for e — > 0, 



(Si-fla)(2-e) = ^-r 2 + 0(e), 

(5i-52)(e) = ^ + 0(e). 

By choosing e sumciently small, we can then guarantee that (51— g2){2 — e) < and (<7i — <?2)(e) > 
0. Hence by deforming the contours if necessary, we may assume that Re (g\ — 52XC) < on 
f 2 U f 8 , while Re (31 — 32) (C) > on ^3 U F7, which ensures that the jump matrices on these 
contours uniformly tend to the identity matrix, exponentially fast as A — > +00. 



3.6 Construction of global parametrix 

Away from the points 2 and 0, we expect that C should be well approximated by the solution 
of the following RH problem, which is obtained from the RH problem l3.8l for C by removing 
all exponentially decaying entries in the jump matrices: 

RH problem 3.9. 

(1) C(°°)(C) is analytic for ( £ C\ ((-00, 0] U [2, 00)). 

(2) C(°°)(C) has the jumps 



cr^C) = c ( _°°\c) 



c^°\o = c [ _°°>(0 



(00) , 



(° 





1 







1 








-1 



















V 


f 1 








°\ 











1 








1 





V° 


-1 





0/ 



for C > 2, 



for C < 0. 



(3) As £ — > oo, we /lave 

c*»> (o = (/ + ^ + o (i) ) dia g ((-cr 1/4 , c 1/4 , (-o 1/4 , c 1/4 ) 



V2 



/ 1 -i 0\ 

1 i 

-i 1 

\0 i 1/ 



(3.30) 



It is worthwhile to point out that C^ 00 ) is independent of A. The RH problem for C 1 ' 00 - 1 can 
be solved explicitly, and its solution is given by 



c(~> (o - diag ((2 - cr 1/4 , c 1/4 , (2 - o 1/4 , c 1/4 ) -J= 



/ 1 -i 0\ 

1 i 

-i 1 

\0 j 1/ 



(3.31) 



where we take the branch cuts of C 1 ^ 4 and (2 — C) 1 ^ 4 along (—00, 0] and [2, 00), respectively. 
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The residue matrix C^ 00 ' in (|3.30[) is diagonal: 



= diag 



(3.32) 



3.7 Construction of local parametrices 

The global parametrix C 1 -- 00 ^ is a good approximation for C only for £ bounded away from the 
points and 2. Here we will construct the parametrices and near these points. Since 
the local parametrices around and 2 can be built in a similar manner, we will only consider 
the local parametrix near 0. Let D(0 : S) be a fixed disk centered at with radius 5 < e, 
and let dD(0, 8) denote its boundary. We look for a 4 x 4 matrix-valued function defined 
in D(0, S) which satisfies the following. 

RH problem 3.10. 

(1) C(°)(C) is analytic for ( e £>(0, <S) \ (R U T 4 U r 6 ). 

(2) For C € D(0, S) n (R U f 4 U f 6 ), has tfie jumps 

Cf(C)=Ci 0) (C)Jc(C), 
where Jc is the jump matrix in the RH problem for C . 
As A — > +oo, we have 



(3.33) 



(3) 



C(°)(C) = C (c 



"4 



(3.34) 



uniformly for ( E dD(0, 5) \ (R U T 4 U T 6 ). 

The solution of the above RH problem for C'- ' can be built via Airy functions and their 
derivatives in a standard way, we follow the theme in [TH1 E] . 
Let j/O; 2/i and 2/2 be the functions defined by 

2/ (C)=Ai(0, y 1 (C)=wAi(wC), lfc(0 = ^ 2 Ai(w 2 0, (3.35) 

e 2 ^ 1 / 3 . Consider the following 4x4 matrix- valued 



where Ai is the usual Airy function and uj 
function ty: 



f 1 




V° 
f 1 




V> 




Vo 








2/o(C) 


2/o(0 


-2/i (0 


-2/i(C) 


-2/ 2 (C) 


-2/ 2 (0 



2/o (0 
1 

2/o(0 



\ 

2/2(0 



2/2(0/ 

\ 

-2/2(0 


"2/2(0/ 

\ 

2/i (0 


2/1(0/ 

\ 

2/i (0 


2/i(0/ 



argCG (0,2tt/3), 



argCG (27r/3,7r), 



argCe (-7T, -2tt/3), 



argCG (-2vr/3,0). 
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The parametrix C^ ) is then built in the following form 

C<°>(0 = S(0* (AC) diag(l, e A92(c) , 1, e~ Aff2(c) ), 



(3.36) 



where E(£) is analytic in -D(0, 5). It is straightforward (cf. [T|)J Sec. 7.6]) to verify that 
given m ([336)1 satisfies items (1) and (2) of the RH problem for C {0 \ To deter mine the analytic 
prefactor E, one needs to use the matching condition (3) and the asymptotics of Airy function 
Ai(C) as C — > oo. A direct calculation gives 



E(0 = C(°°\Q 



(I 














1 

—iy/n J 



(AC) 1/4 


Vo o 



\ 







1 

o (ac)- v v 



This completes the construction of the local parametrix . 

3.8 Final transformation: C >-)■ D 
Definition 3.11. We define the final transformation 

C(C)(C (0) ) -1 (C)> in a ^-neighborhood of 0, 
D(C) = <( C(C)(C (2) ) -1 (C), in a ^-neighborhood of 2, 
C(C)(C (oo) ) -1 (0, elsewhere. 



(3.37) 



From the construction of the parametrices, it follows that D satisfies the following RH prob- 
lem. 

RH problem 3.12. 

(1) D is analytic inC \ Y,jj, where Ed is shown in Figure^ 

(2) D has jumps -D+(C) = D_(C)Jd(C) f or C G ^d, where 




C^\0(C^)-H0, on\{\=5, 
CM(C)(<70>)-i(C), on|C-2|=5. 
C (oo) (C)^c(C)(C (oo) ) _1 (C), o« tiic rest o/E D . 



£>(C) = /- 



C 



c 2 



(3.38) 



(3.39) 



The jump matrix for D satisfies 

J D (0 = I + O(l/X), 



as A — > +oo, 



uniformly for z on the circles |z| = <5 and |z — 2| = 6, and the jumps on the remaining contours 
of E^ are exponentially close to the identity matrix. In particular, we note that 

MO =I + e A ^-^)(0 C (°°)(C) (E 3A - E 2 .j) (C^)" 1 ^), for C G f 2 U f 8 , (3.40) 
MO =I + e-^-^C^HC) (Ei,2 - ^4,3) (C^yHO, for C G f 3 U f 7 , (3.41) 



25 



Figure 4: Contour Yjo for the RH problem 13. 121 for D. 



and 



-2Aai(0 





(2-C) 1/2 

\ o o 



o (2-cr 1/2 o\ 



-i 

0/ 



-2A 92 (C) 



(0 








— i 








V» 


£l/2 



-A(si+ff 2 )(0 







(¥) 



1/4 



^- (C(2-0) 



1/4 



-*fe) 1/4 (C(2-C))- 1/A 

o (c(2-or 1/4 o 

- (C(2 - c)) 1/4 o 

i(^) i. / 



1/4 



2-C 

c 





1/4 



for C S (e,2 - e), with the aid of ([3351) and (l3~3Tj) . 

By standard arguments in jTBl HZ] , we then conclude that 



D(()=I + 



A(ICI + i) 



as A — > +oo, uniformly for (gC\Sd. 

It follows from f|3 . 27[) and the above constructions that 

c = -Vsi(-E>i)i,3 + s?, c = -i>/si(-Di)2,4) = -i-\/si(- D i)i,4 ! 



(3.42) 



(3.43) 



(3.44) 



where we used (|3.37[) for large £ and the fact that C^°°' ) in (I3.32[) is diagonal. 



3.9 Proof of Proposition 13.11 
Solvability of the RH problem for si large 

In the above, we applied a series of invertible transformations M \-> A>- > Bh ;> C H> D, so that 
the matrix-valued function D exists and uniformly tends to the identity matrix as s\ — > +oo. 
This immediately implies the solvability of the RH problem for M for si sufficiently large. 



2G 



Asymptotics of c, c and d as si — > +00 



3/2 



it follows from (pQgj) that Di = 0(A _1 ) = 0(5! 3/2 ) as si -> +00. Then by 



5 ? + o( s r 1 ), c^o^r 1 ), 



o^r 1 ), 



Recalling A 
fl33I|), we have 



as si — >• +00. This, together with our assumption (|3.4p . implies p.2[) and Q3.3p . 

It remains to establish (|3.ip . This requires more effort. First we make some observations on 
the jump matrix Jd in (I3.38[) . The matrix Ju has the 'checkerboard' sparsity pattern 



/* * 0\ 

* * 

* * 

\0 * */ 



(3.45) 



on the circles \z\ = 6 and \z — 2| = 5, on (—00, e) U (2 — e, 2) and on the contours f 1 uf 9 Ul^ UIV 
Here, * denotes an unspecified matrix entry. 

Let be the solution to the following RH problem. 

RH problem 3.13. 

(1) is analytic in C\ £d, where is shown in Figure^ 



e -2A 31 (0 



„-2Ag 2 (C) 



-HQJdMC) for(€X D , wh 


ere 


( i (2- 


C) -l/2 


°^ 











(2-C) 1/2 


—i 





\ 







/o \ 






Q -i C~ 1/2 









3 




\o c 1/2 » y 







I, 

MO, 



/orCe(e,2-e), 



/or cer 2 ur 3 ur 7 ur 8 , 

elsewhere. 



(3) ^4s £ — > 00, we ftare 



(i) 



< 2 



(3.46) 



It is easily seen that J D (\) still tends to the identity matrix as A — > +00, hence, there is a 
unique solution ZjW to this RH problem if A is large enough. Moreover, we have 



£>W(C)=I + 



A(|CI + 1) 



(3.47) 



as A — > +00, uniformly for (gC\Sd. 

Since all the jump matrices for D^> have the checkerboard pattern (|3.45[) . the solution D^> 
has the same pattern. Obviously, the residue matrix in (|3.46l) also has the pattern (I3.45p . 



Definition 3.14. We define 



r(0 = d(0(dW)- 1 ((). 



(3.48) 
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Figure 5: Contour for the RH problem 13. 151 for R. 



The matrix-valued function i?(C) satisfies the following RH problem. 
RH problem 3.15. 

(1) R is analytic in C \ S/j 7 where S/j is as in Figure [5j 

(2) R has the jumps R+(C) — R-(C)Jr(C) f or C € where 

MO = D^iOMOi^rHO- (3-49) 

(3) ^4s C — >• oo, we have 

R(C)=I+^+o(±y (3.50) 

Since R = Dpi 1 ))" 1 , we obtain from (pTMj) . (gUSjl and (|3"30)) that 

R l =D l - D<P. 

This, together with (13.441) and the checkerboard pattern of D± , gives us 

d= -iViT(i?i)i,4. (3.51) 

Next, we give some estimates for Jr defined in (|3.49p on By (|3.40[) . (|3.41[) and the above 
definitions, it is readily seen that 

Jr(0 =1 + (V (91 ~ 92)(C) ) , uniformly for ( G f 2 U f 8 , 
Jr(0 =1 + f e -A(fll-fl2)(c) ) , uniformly for ( € f 3 U f 7 , 

as A — >■ +oo. 



2s 



On the interval (e, 2 — e), we have 

MO = D^iOMOJnlM^rHQ 

with 

o -*(^c) 1/4 

(¥)'" 



-A( gi + S2 )(C) 



1/4 



\- (C(2-C)) 1/4 

x (J+ (9(e- 2Asi (^) + 0(e- 2X92 ^)), 
as A — > +oo; see (|3.42p and the definition of Jjju). Hence 

e -A(gi+g 2 )(0 



(C(2-C)) 
o (C(2-c))" 1/4 o 
(C(2-c)) 1/4 o 





2-C 



1/4 



1/4 



Jfl(C) = -f + 
/ 



-^ a) (C) 



o 



1/4 



1/4 



^(2^)' o (C(2-C)) 

o (C(2-c))" 1/4 o 





'(¥) 

o -(C(2-C)) 1/4 
(C(2-C)) 1/4 o 

x (/ + 0( e - 2Asi ^)) + 0(e- 2As2 ( c )))(D ( L ) )- 1 (C) 



1/4 



«(¥) 





-i/A 



1/4 



for (e (e,2-£). In particular, it then follows from p. 471) and the above formula that, for the 
(1,4) entry, 

e -A( S i+ S2 )(C) 



(^)i, 4 (C) = 



(l + OtA- 1 )), (€(e,2-e), 



(3.52) 



2(C(2-C)) 1 ¥ 
as A — > +oo. 

Now we observe that the function g\ + g2 has a unique minimum on the interval (0, 2), which 
is attained at the point 

C*:=^— G(0,2). (3.53) 

To "T 1 



Indeed, straightforward calculations using (|3.20[) and (I3.21[) yield 

(.9i +52)(C ) 



3 + 1 



(3i+.92)'(C) = 0, 



(.9i+.92)"(C) 



v/2 (r 2 + l) 3 /2 



4 r 2 

We may assume that the number e > is small enough so that 

(.91 + .92)(C) < Re (g 2 - 9l )(C), for C € f 2 U f 8 , 



(3.54) 
(3.55) 
(3.56) 



29 



and 

(3i + 92)(C) < Re ( 91 - g 2 )(0, for C € f 3 U f r . 

Thus the exponents in the exponential estimates (13.40^ — f|3.41|) are strictly smaller than — X(gi + 
92) (C*) with (|3.54l) . Therefore, the jump matrix Jr behaves as 

MO =/ + 0(e- A(si+S2)(r) ) , A -» +00, (3.57) 

uniformly on Sr. From p.57[) . standard theory implies a similar estimate for R itself: 

i?(C) =/ + 0(e- A(9l+92)(c * ) ) , A — > +00, (3.58) 

uniformly for £ € C \ Sr. 

Since i? + = i?_ + (Jr — J) + — I)(Jr — I), it then follows from the Sokhotski-Plemelj 
formula that 

m=i+-f {jR - I] } x) ^-l v^WtzlMte 

= I + — I ^— ii ^.da; + O , A -> +00, 



27r * J-£ R x-C V c 

on account of (|3.57j) and (|3.58[) . For the (1,4) entry, we get 

{Rx) lA = -^- i J^ (J R ) 1A (x)dx + 0[e- 2X ^ +9 ^ c ^, (3.59) 

and by ([3~5"Tj) . 

( Jji)i,4(x)da? + O (V 5 i e ~ 2A( * 1+fla)(n ) , A — > +00. (3.60) 



2tt 



The main contribution to the integral in (|3.60p comes from a neighborhood of £*. For any 
given S > we have 

as A — > -f-00 , by virtue of (|3.52j) . Now a standard saddle point approximation (Laplace method, 
cf. g0l|46]) yields 



/IT / 27T p -A(9l+ff 2 )(C*) 



as A — > +00. Inserting p.53p . (|3.54[) , p. 561) and A = s^ 2 into p.6ip . we finally obtain 

2,1 2 , 1 \ V 4 4V2 r 2 .3/2 

i±ii±^ e 3 (l + O^ 2 )) 



V2 f 4 r 2 A 1/2 /rf+lrf^ 174 ^ — ^ 
40F S } /4 ' l^(r 2 2 + l) 3 / 2 J ' V 2 2r 2 2 



20Fsl /4 V2(r 2 2 + 1 



1/4 4n/2 r 2 .3/2 

3 /„2 1 1 1 , . ; ■> . 

Si — » +00, 

which is p. II) with r\ = 1 and s 2 = 0, as desired. 

This completes the proof of Proposition 13. II □ 
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4 Proof of Theorem Q 



In this section we prove Theorem 12.31 on the solvability of the RH problem for M (£) by using 
the technique of a vanishing lemma. To this end we follow basically the scheme laid out in 
[T71 H7] , although the argument is somewhat more involved because our RH problem is of 
size 4x4 whereas the usual dimensions treated in the literature is 2 x 2. 

Following [17], the proof consists of three steps, which as in [T7J are called Step 1, Step 2 and 
Step 3. 

Step 1: Fredholm property 

Standard theory show that the RH problem for M is associated to a singular integral operator. 
The first step is to show that this operator is Fredholm. To this end we are going to apply a 
series of transformations M *-> M^> i— > M^ 2 \ The first transformation M \— > AfW is defined by 

M«(C) = M(C)e(C), (4.1) 

where 

9(C) = diag (>«, e e ^\ e~ e ^\ e~ 9 ^) . (4.2) 

This transformation will kill the exponential factor in the asymptotics (|2.37[) . at the expense of 
complicating the jump matrices. 

The second transformation \— > is defined by 

M^(C; r u r 2 , Sl ,S2) = M (1) (C; n, r 2 , s u s 2 )(M«)- 1 (C; r 1 ,r 2 ,s* 1 ,s 2 ), (4.3) 

where is a fixed but sufficiently large real number for which we already know that M and 
therefore exists; see Proposition 13. II Then M@' satisfies the following RH problem. 

RH problem 4.1. 

(1) MW(C) is analytic for f e C \ (uLi r * U (jLe r *) ■ 

(2) For C e r fe; M< 2 !(C) has the jump 

m| 2) (0 = m1 2) (C)4 2) (0, 

where 

4 2 \() =Af W (C;r 1 ,r 2 , S I, S2 )e: 1 (C)^e + (C)(M| 1) )- 1 (C;r 1 ,r 2 ,^, S2 ), 
and where Jk is the original jump matrix on Tk denoted in Figure 

(3) As C oo, we We M< 2 )(C) = / + 0(C _1 ). 

Since the transformations M h-> i-> M^ 2 ) are invertible, we have that the original RH 

problem for M is solvable if and only if the one for is solvable. 

(21 

One checks that the jump matrices Ji (£) converge to the identity matrix whenever £ — > oo 
or £ — >■ 0, along any of the 8 rays 1^, k = 1, . . . , 4, 6, . . . , 9. In other words, the jump matrices are 
normalized at infinity and continuous at zero. Then it follows from the techniques of Deift et al. 
[17] that the singular integral operator associated to the RH problem for M^ 2 \C,) is Fredholm. 
The same statement then holds for the original RH problem for M(Q. 



31 



Step 2: the Fredholm index is zero 

The next step is to show that the Fredholm index of the Fredholm operator associated to M^> is 
zero. This follows by a continuity argument in the same way as in |17j . Indeed, we have that the 
Fredholm index is a continuous function of s± which takes only integer values. For large values 
of si, the Fredholm index is equal to zero, and therefore this must hold for all s±. 



Step 3: triviality of solution to the homogeneous version of the Riemann-Hilbert 
problem 

The third and final step is to show that the 'homogeneous' version of the RH problem has only 
the trivial solution. This is also known as a vanishing lemma. As in 17, Page 1402], this reduces 
to studying the following RH problem, which is the homogeneous version of the RH problem for 
M(C). 

RH problem 4.2. We look for a 4 x 4 matrix-valued function Af( 3 )(£) satisfying 

(1) M( 3 '(C) is analytic for ( G C \ (uLo r k 

(2) For C G T k , M^(C) has the same jumps as M{Q, see Figured 

(3) As £ — > co, we have 

M( 3 )(c) = (otc 1 )) dia g ((-cr 1/4 , c 1/4 , (-o 1/4 , c 1/4 ) 



/ 1 -i 0\ 

1 1 i 

V2 ~i 1 

\o i o iy 



e- 4 (C), 



(4.4) 



where we recall the notations (|4.2I) and fl2.38|) — (|2.39|> . 

(4) M< 3 )(C) is bounded near ( = 0. 

Note that the only difference between the RH problems for M(£) and M( 3 )(£) is in the 
leftmost factor of the asymptotics; compare (|2.37[) with (|4.4[) . 

We need to show that the only solution to RH problem 14.21 is when M^> is the zero matrix. 

First we apply a contour deformation to bring all jumps to the real axis. Denote by f2/c, 
ft = 0, . . . , 9 the region between the rays Tk and Tk+i in Figure [2J with Tio := r . We define 

MW(C) forCeC\ULo r fc b y 



M( 4 ) = 



M< 3 


Ji J2, 


for C G O 


M( 3 




for (efli 


M< 3 




for (en 2 


M< 3 




for (efi 3 


M< 3 


T-1 T-1 

J 4 ^3 ' 


for C G 4 


M< 3 




for (efi 5 


M< 3 




for (efi 6 


M< 3 


1 


for (e 7 


M< 3 




for (efi 8 


M< 3 


T-1 T-1 

J 9 J 8 > 


for (efi 9 



(4.5) 
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where Jk, k = 0, . . . , 9, denotes the jump matrix for M(£) on I\. in Figure [3J Next we set 

Af( 4 )(C) = MW(C)6(C). 

RH problem 4.3. The 4x4 matrix-valued function M^ 4 \Q satisfies the following RH problem: 

(1) Af (4) (C) is analytic for ( 6 C\M. 

(2) For x € K, we /lave the jump relation 



(I 





1 


A 





1 


1 


1 








1 





V» 








V 



M| 4) (x) =M l i ) (x)ez 1 (x) 



where the orientation of the real axis is from left to right. 
(3) As ( —} oo, we have 



e+(x), 



(4.6) 



(4.7) 



M (4) (c)-o(r 3/4 ). 

(4) Af( 4 '(C) is bounded near ( = 0. 
Proof. The fact that M W(C) does not have jumps on Tfe, k = 1, 4, 6, 9, follows immediately 



from (|4.5p . The jump matrix on the real axis (|4.6p follows from the relations 
Je,JgJoJ\J2 — (J3J4J5J6J7) 1 = 



fl 1 l\ 

111 

10 

\0 1/ 



also note that we reverse the orientation of the negative real axis. As for the asymptotics (|4.7p . 
one should check that 6 _1 J\ J2O is bounded for ( G fi , Q~ 1 J 2 & is bounded for £ e fix, and so 
on. Equivalently, e® 1 ^ should be bounded for £ 6 f2o, e 01 ^' should be bounded for £ € fii, 
and so on. These statements are easily checked from (|2.38[) - (|2.39p and (|2 . 33|) . 

Next, we define a new 4x4 matrix- valued function M( 5 )(£) by 

M( 4 )(£), C m lower half plane of C, 

M(4) (C) f t _ rf 2 ) > C in upper half plane of C. 



M< 5 )(C) = 
Then M^(() has the jump 



□ 



(4.8) 



Mf(x) = M ( _ 5) (x)J(x), 



x e 



(4.9) 



with 



J(x) := ez\x) 



ez\x) 



fl 





1 
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1 










1 





e 
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-1 


°\ 
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1 







-1 
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V° 
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0/ 




h 



ez H (x), 



(4.10) 
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and where the superscript H denotes the inverse Hermitian conjugate. Similarly, we will use H 
to denote the Hermitian conjugate. 

Now define a new 4x4 matrix-valued function Q(Q = M (5 ) (C)(M( 5 ) {()) H . Then Q(() is 
analytic in the upper half plane of C and it decays with a power £~ 3 / 2 as £ — > oo. A standard 
argument based on contour deformation and Cauchy's theorem shows that 



Q+(x) dx = 0. 

JS. 

Hence, 

/ R M^ W( M^))^ = 0. 
By adding this relation to its Hermitian conjugate, we find 

Mi 5) (x) (J(x) + J H {x)) (M { _ 5 \x)) H dx = 0. 

But from (|4.10l) we have that 



J(x) + J H (x) = 29 



= 2 



/l 1 0\ 

110 



\0 0/ 

e -0i,-(V)\ 









-H 



V o / 

Substituting this expression in (|4.11[) yields 



M { _ 5) (x) 





V o ) 



L-ei,-[*) g-fti.-W o o) (Mi 5) (a;)) H ds = 0, 



which obviously implies that 



M ( _ 5) (x) 



/ e -*i,-(*)\ 

-03,-0) 


V o / 



0, a; e 



(4.11) 



(4.12) 



(4.13) 



Inserting (|4.13p into the RH problem for M^ 5 \ we see that the jump relation (|4.9[) - (|4.10j) reduces 
to 








-1 


°\ 











-1 


1 













1 





0/ 



m| 5) (x) = M (5) (x)ez 1 (x) 



By tracing back the transformation *-> M^ b \ we find 

m^\x) = M {i) (x)ez 1 (x)e+(x). 



ez H (x). 
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This now decouples into four scalar RH problems for the individual column vectors. For each of 
these scalar RH problems, one can use an argument based on Carlson's theorem jTTJ Page 1406] 
to conclude that it has only the zero solution. This then implies that also must be the zero 
matrix. By tracing back the transformation \— > the same conclusion holds for 

Thus we see that the homogeneous version of the RH problem for M indeed has only the trivial 
solution. This ends the proof of Step 3, and therefore the proof of existence of M by means of 
the vanishing lemma. □ 



5 Proof of Theorem 12.41 

In this section we prove Theorem 12.41 on the Painleve II behavior of the numbers c, c and d in 
5.1 Symmetry properties 

We will need a few properties of the RH problem 12.21 The first property follows from exploiting 
the symmetry of the problem. In what follows, we use the elementary permutation matrix 



J 



1 

1 



and recall that 1% denotes the 2x2 identity matrix. 

Lemma 5.1. (a) For any fixed r\, r 2 > and s±, s 2 G we have the symmetry relations 

"®-{o -/>«>( o -*)■ (5 ' 1) 

where the bar denotes complex conjugation and 

M - T «> = (_" 2 o) M(0 U ~°)' (5 ' 2) 

where the superscript ~ T denotes the inverse transpose. 

(b) Denoting with M(£; ri, r%, s\, S2) the solution M(£) with parameters r\, T2, s\, S2, then we 
have the relation 

Af(-C;ri,r 2 ,si,s 2 ) = °j^J M((;r 2 ,r 1 ,s 2 ,s 1 ) . (5.3) 

Proof. One easily checks that the left and right hand sides of (|5.ip satisfy the same RH problem. 
Then (|5.1[) follows from the uniqueness of the solution to this RH problem. The same argument 
applies to (pT2"j) and (jO)) . □ 

Corollary 5.2. For any fixed r\, r 2 > and s\, s% S M., we have 

MS XX -I 



and 
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Consequently, Mi takes the form 



Mi = 



fa b ic id\ 

-b —a id ic 

ie if —a b 

\if ie —b a/ 



where a, a, b, b, c, c, . 
We also have 



are real constants depending parametrically on ri, r 2 , si, s 2 . 



-Mi(ri,r 2 ,si,s 2 ) 
a nil consequently we have in (|5.4p that 



a(ri,r 2 ,s 1 ,s 2 ) = a(r 2 , r x , s 2 , Si), 
c(ri,r 2 ,si,s 2 ) = c(r 2 ,ri,s 2 ,si), 

Proof. This follows from (|2.37[) and Lemma 15.1 



Mi(r 2 ,ri,s 2 ,si) 



6(n,^2,si,s 2 ) = 6(r 2 ,ri,s 2 ,si), 
e(ri,r 2 ,si,s 2 ) = e(r 2 , n, s 2 , Si). 



(5.4) 



(5.5) 
□ 



5.2 Lax pair equations 

We next obtain linear differential equations for M with respect to both £ and si, s 2 . This system 
of differential equations has a Lax pair form and the compatibility condition of the Lax pair will 
then lead to the Painleve II equation in Theorem 12.41 
We start with a differential equation with respect to (. 



Proposition 5.3. We have the 



erential equation 
dM 



UM, 



(5.6) 



U := 



[Mi,iriE 3t i - ir 2 E iy2 ] 



(5.7) 



where 

( in 0^ 

ir 2 

i(rxC-Si) 

\ ~i(r 2 C + s 2 ) / 

and where [•, •] denotes the commutator. 

Proof. Since the jump matrices in the RH problem for M do not depend on we have that 
has the same jump properties as M. It follows that U :— ^M-M^ 1 is entire. From the asymptotic 



behavior of M in (|2.37[) . we have as £ — > 00, 
1 

4C 

/ 




u 



— diag(-l, -1, 1, 1) + ( I + + 



Mi 



M 2 



i(ri - siC x ) 


i(nC-si) 

V -i(r 2 ( + s 2 ) 

Mi Ml - M 2 





«(r 2 + s 2 C- 1 ) 





x II- 



c 



( 2 



+ o(c z )- 



(5.8) 



Thus the entries of U are polynomial in ^ with degree at most one. Dropping the non-polynomial 
terms in (|5.8p . we obtain (15.71) . □ 
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The proof of Proposition 15.31 also yields the following. 
Lemma 5.4. The entries a,a,b,b, . . . of Mi in (|5.4[) satisfy the identities 

(2a + c 2 )n = r 2 d 2 + s u (5.9) 
(2a + c 2 )r 2 = rid 2 + s 2 , (5.10) 
r 2 b — r\b = (r 2 lz — r\c)d. (5-11) 

is equal to zero. For the upper right 2x2 block, we then 



Proof. The coefficient of £ 1 in ((5 
obtain by using (|5.4[) that 



-is\ 
is2 



a b 
—b —a 



iri 






ir-2 



iri 
ir 2 



—a b 

—b a 

c d 

d c 









-ir 2 



c d 
d c 



These are four identities and they give us (|5.9[) - (|5.11[) . 



□ 



In the special case r\ — r 2 and s\ — s 2 , we see from (|5.5[) that a = a, b = b, and so on. The 
equation (15.11[) then reduces to = 0, while (|5.9p and (|5.10p are the same. So in that case, the 
system (|5.9[) - (|5.11[) reduces to the single relation 



(2a + c 2 )r = r 2 d 2 + s, 



with r = r± = r 2l s = s± = s 2 . 



(5.12) 



We obtain more differential equations by taking a derivative of M with respect to the param- 
eters Si, S2- 



Proposition 5.5. We have the differential equation 

dM 



dt 



VjM, j = 1,2, 



where 



Vi = 



Vo = 



( 





-2i 
















— 2«C 











V o 








0/ 


/o 








°\ 











2i 
















2< 





0/ 



[Mi, -2^3,1 



(5.13) 



(5.14) 



[Afi,-2iB 4 , 2 ] 



(5.15) 



Proof. Since the jumps in the RH problem for M do not depend on sj, we have that dM/dsj 
has the same jumps as M and so Vj — (dM /dsj) M _1 is entire. As ( — > oo, we find for j = 1, 



Vi = [I 



Mi 

T 



/ 



-2iC 

\ o o 



-2i 0^ 





0/ 



7- 



Mi 



(5.16) 



Keeping only the polynomial terms in £, we obtain (|5.14l) . The other relation f|5 . 1 5[) is proved 
similarly. □ 
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5.3 Compatibility and proofs of (|Z^I]) - (|Z33j) 

The compatibility condition for the two differential equations (|5.6p and (|5.13p is the zero curva- 
ture relation 

^- d ^ + U Vj - Vj U = 0, j = l,2. (5.17) 
This is, in view of ([BTf]). (j5T4|) and (j5T5l) . 

A [Mi,iri£7 8 ,i - ^2^4.2] = -iE 3 ,i + , (5.18) 

osi 

^- [Mi,iriE 3 i - ir 2 E ia ] = -iE ia + [V 2 , U] . (5.19) 

OS 2 

Both (|5.18p and (|5.19[) give us a system of 16 differential equations for the entries of M\. 

With the above preparations, we are ready to prove (I2.4ip ~ (|2.43p in Theorem 12.41 concerning 
the Painleve II behavior of the numbers c, c and d. To this end, we work with (|5.18l) and first 
derive the differential equation satisfied by d. 

The entries in the matrix relation (|5 . 18[) can be obtained from a lengthy calculation or with 
the help of a symbolic software package such as Maple. For the (1, 1) and (2, 2) entries of (|5.18p . 
this yields 

ri |i = 2r 2 d 2 + 2 Sl , ^- = 2d\ (5.20) 

OSi OSi 

which imply 

d(r lC ~ r 2 c) 

5 =2si. (5-21) 

The (1, 2) and (2, 1) entries of (|5.18p give expressions for the partial derivative of d: 

3d — To 

— = 2cd- 2b = — (Zed- 2b), (5.22) 



and the (1,4) entry gives 

d(rib + r 2 b) 



r 2 (4ad + Aad + 2bc + 2bc) + 2s 2 d. (5.23) 



ds i 

On the other hand, the identity (|5.1ip implies 

ri(r\b + r 2 b) = (r\ + r 2 )b + r 2 {r Y c - r 2 c)d, 

which upon differentiation leads to 

d(rib + r 2 b) . 2 2 . db , „ dd , ntS 

ri^^ = (r? + r\ — + r 2 {r x c - r 2 c) — + 2r 2Sl d, 5.24 

OSi OSi OSi 

where we have made use of (|5.2ip . Combining this with (|5.23l) yields 

db dd ~ 
( r i + r i)^ 1" r 2 {r x c -r 2 c\- h 2(r lS2 + r 2 s 1 )d = -r ir2 {Aad + 4ad + 26c + 26c). (5.25) 

We next eliminate a and a from the right-hand side of f|5.25[) with the help of (|5.9p and (|5.10p . 
This gives us 

/ 2 2\ 9b . „ dd „, . , 

( r i + r 2)-^ h r 2{ric - r 2 c)- h 2(ris 2 + r 2 si)d 

osi osi 

= -2(rf+rl)d 3 - 2(r 1 s 2 + r 2 s 1 )d + 2r 1 r 2 (c 2 d-bc + c JZ d-bc). (5.26) 
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We move the last term in the left-hand side of (|5.26p to the right, and rewrite the last term in 
terms of dd/dsi by using (15.221) and (|5 . 1 1[) . It then follows that 

. 2 2\ 9b . dd 
( r i + r 2)x \-r2(riC- r 2 c) — 

= -4(ns 2 + r 2 s 1 )d - 2(r 2 + r 2 2 )d 3 + n(r 2 c + r{c) — , 

OSi 



or equivalently, 



- c# = - 4(ri5 2 2+r 2 2Sl) d - 2d\ (5.27) 



Taking a derivative of the second identity in (|5.22[) with respect to si and using (|5.27[) to 
eliminate db/dsi — cdd/ds\, we obtain 

d 2 d r 2 ( dc 4(ris 2 +r 2 si) 



2- 2 ^[^- d + J^J J d+2d*). (5.28) 



ds\ ri \ds\ r\ 



2 



This, together with the fact that dc/dsi = 2d 2 (see (|5.20p ). implies 

^ = 2^ f2d 3 + 4(riS 2 2 + r 2 2Sl) rf + 2d 3 
asf n \ rj + r 2 

= 8 r 2 (r lS2 + r 2Sl ) d + 8 r 2d 3 

r"i r;[ + r§ ri 

The differential equation (15.29[) is a scaled and shifted version of the Painleve II equation. 
Indeed, we have that q satisfies q" = sq + 2q 3 , if and only if 

f(s) = aq(f3s + j) 

satisfies 

f" = ^(/3s + 1 )f + 2^f. (5.30) 
Comparing this with (|5.29[) . we see that we need a, (3 and 7 so that 

r\ a z r\ r\ + r 2 r% r\ + r 2 



which means 



( M/6 o 2/3 2/3 

_ (^ir 2 ) 1/D _ 2r 2 ' 2r x * s 2 

a ~ /„9 , JM /?! ' P 



(r 2 +r2)V3' " rJ /3 (r 2 +r2)V3' r V 3 ( r 2 + r 2)i/ 3 ' 

Therefore, we have proved that 

d = aq ( Psi + 7 ) = (r2+r2)1/3 g ( v (rir2) l/3 (r 2 +r 2 )1 /3 J ( 5 - 31 ) 

with g being a solution of the Painleve II equation, which is (|2.41[) . The fact that q is the 

Hastings-McLeod solution follows from the asymptotic behavior of d in (|3.ip . To see this, we 
rewrite Q3.1[) in the following form: 

(rir 2 ) 1/6 ../ 2(ris 2 +r 2 si) \ _ 3/2 

rf = /„2 I ^2M/3 Al /I I M/3/-2 I _2M/3 ( X + ( S 1 ))> S l "> +0 °' ( 5 " 32 ) 



(rf + rf) 1 ^ y ( rir2 )i/3( r 2 + r 2)i/ 3 
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by virtue of the asymptotics of Airy function; cf. [T]. Comparing this expressions with (|5.31l) . 
we see that q(s) ~ Ai(s) as s —> +00, and so q is the Hastings-McLeod solution. 

To establish (|2.42[) and (|2.43l) . we first derive an identity for the si-derivative of c and c, 
respectively. Noting that u' — —q 2 , where u is the Hamiltonian in (|2.44l) . the first equation in 
(|5.20[) can be written as 



dc 



2^d 2 



9 r 2 (rira) 1 / 3 2 
n (r 2 +r|) 2 /3 9 



d 

dsi 



2/3 

f 1 



2(r 2 si+ri5 2 ) \ 
A /3 rl /3 (rf+rl)^J 

( 2(r 2 s 1 +r 1 s 2 ) 



''1 



ri /3 (rf + r|) 1 /3 ^v^' 8 ^ + rf^y ri 



1/3 1/3, 



Similarly, it follows from the second equation in (|5.20[) that 



= 2d 2 = — 



2/3 



2(r 2 si 



7"lS 2 j 



(5.33) 



(5.34) 



By integrating (|5.33p and (I5.34[) with respect to sx, respectively, we then obtain (12.421) and (|2.43[) 
on account of (|3 ,2|) — (|3 ,3[) and the fact that u(s) = o(l) as s —> +00 if q is the Hastings-McLeod 
solution of the Painleve II equation. 

This completes the proof of Theorem 12.41 □ 



Part II 

Non-intersecting Brownian motions at 
the tacnode 

Remark on notation and conventions: Throughout the next sections we work under the 
assumption of the double scaling limit as described in Section 12.31 Throughout most of Sections 
E] and [7] n will be large but fixed. Most notions depend on n, such as pj, a,j, bj, ctj, /3j, Vj, 
Xj, and so on, although this is not indicated in the notation. Their limit values as n — > 00 are 
denoted with a star, such as p*,a*,b*. Also recall that the temperature T = 1 and that the time 
t = £ C rit is fixed (independent of n) according to (|2.27[) . 

By the translational symmetry of the problem, it will be sufficient to give the proofs in case 
where the tacnode is at the origin. That is, we will assume that 

£crit := P2 — a * — °- 
From ([2Tn?|) - (|2~20")) . this yields the relations 

since T = 1. For ease of notation we also set 

a :=-a 2 , (f := p x . 



(1 - t)a* 2 - tb* 2 = 2yfp* 2 t(l-t), (5.35) 
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Then the Brownian paths are asymptotically distributed on the two touching intervals [— a*,0] 
and [0,0*] where (use (jQ5j) to simplify (|2~I%j) - (12^2111) ') 



a* :=4Jp*t{l-t) > 0, (5.36) 



P* :=^p\t(l-t) > 0. (5.37) 
We will use these conventions throughout Sections [SHE1 

6 Modified equilibrium problem, ^-functions and A-functions 
6.1 Modified equilibrium problem 

In the steepest descent analysis of the RH problem 12.11 for Y, we will use functions that come 
from an equilibrium problem for the two external fields 

Vj(x) = (x 2 - 2((1 - t)aj + tbj)x) , j = 1,2. (6.1) 

In the usual equilibrium problem with external field (16.11) . one asks for a measure that mini- 
mizes the energy functional 

J J log \ x - y \ d K x ) d Kv) + JVj(x)dn(x) 

among all measures on the real line with mass J d/i — pj. Since Vj is quadratic, the solution is 
a semicircle law with density 



2nt( l ~tW ^ ~ ~ l£ [ a J'fti ( 6 - 2 ) 

with 

a.j = (1 - t)a,j + tbj - 2Jp~t{l - t), (6.3) 
Pj = (1 - t) aj + tbj + 2^J Pj t(l - t); (6.4) 

see also ([239]) -([23T ]l . 

The limiting situation corresponds to al,a*,,b*,b*,. Then a* 2 < P* — = a\ < P* and the 
two semicircles (16.21) meet each other in the origin. For finite n however, this may not be the 
case. Indeed the two semicircles may be separated, or may overlap depending on the situation. 
In what follows we modify the equilibrium problems in such a way that the two minimizing 
measures have supports [0,P] and [a, 0], respectively, that meet at 0. We allow the measures to 
become negative, so we will be dealing with signed measures. The modification of an equilibrium 
problem to prepare for a later steepest descent analysis of a RH problem was first done in [IB] , 

We assume that a-x < < pi, which is certainly the case if n is large enough. 

Definition 6.1. In the modified equilibrium problem for V% we ask to minimize 

log- -dp,(x)dfj,(y) + I Vi(x)d/J,(x) 

\x-y\ J 
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among all signed measures y, on [0, oo) with total mass J dfi — pi and such that /i is non-negative 
on [ai, oo), where ct\ is given by (|6.3I) . 

In the modified equilibrium problem for V2 we ask to minimize 

log] rd(j,(x)d(j,(y) + [ V 2 (x)dn(x) 

\x-y\ J 

among all signed measures [i on (— oo,0] with total mass J d[i = P2 and such that /x is non- 
negative on (—oo,^! where 02 is given by 



We can explicitly find the minimizers for these modified equilibrium problems. 

Proposition 6.2. (a) The minimizer in the modified equilibrium problem for V\ is the signed 
measure \i\ on [0,/3] with 



(3 = <*i + ft + V"i + ffi - <*ift > (6 5) 



given by 



dx 27rt(l — t) V x 
with 



dfMl X - &1 x€[O,0\, (6.6) 



ai+Pi- \Jol\ aifh 
01 = g ■ (b-'J 

(b) T/ie minimizer in the modified equilibrium problem for V2 is the signed measure /12 on 
[—a, 0] with 



a = ~"2 - ft + 2y/dj + 0j- a 2 >Q 



given by 



with 



d[i2 5 2 -x x + a 

~d^ = MR V~ * G[ ~ a ' 0] ' (6J) 



3 

Proof. Let #i be the ^-function associated with ^i, i.e., 

Si(z)= [log{z-s)d t x 1 (s), zeC\(-oo,/3]. (6.11) 



Clearly, gi is defined and analytic in C\ (— 00, /3]. It then follows from the variational conditions 
of the equilibrium problem that 

g 1 ,+ (x)+gi,-(x) = Vi(x) + c 1 , x € (0,0), (6.12) 

for some constant ci depending on n. Differentiating both sides of (|6.12p with respect to x gives 

ffL+(aO+si,-(aO = VY(aO, ( 6 - 13 ) 

or equivalently, 

=-Qn'W-5iW) (6-14) 
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for x e (0, 8). This, together with the fact that gi t+ and g± t _ only differ by a constant for x < 0, 
implies that (^V{(z) — g'^z)^ 2 is an entire function in the complex plane. Since this function 
only has polynomial growth for z large and fix is supported on [0,/3], it is readily seen that 

, ., A 2 (z - <5i) 2 (z - 8) 



for some 5± < 8. Hence, 



% = ~Im [^V;(x)-g[(x)) (6.16) 



da; 7r 



+ 



■'' //? — , x€[O,0\, (6.17) 



27rf(l - i) 
which is (|6.6I) . 

To obtain the representations of /? and Si, we expand (|6.15l) as z — > oo. Comparing the 
coefficients of order 0(z 2 ) and O(z) on both sides and taking into account (I6.3p - (|6.4p leads to 

+ + ft, 

*i(2 ) 8 + « 1 ) = ai J 8i. 

Solving the above algebraic equations for <$i and 8 gives us (|6.5p and (|6.7p . 

The explicit formula for /i2 stated in item (6) can be proved in a similar manner as for 
To that end, we need to use 



*(^ 2 '(*)-^(*)Y 



(z-J 2 ) 2 (z + Q ) 



where 

<? 2 (z) = / log(z-s)dn 2 (s), 2€C\(-oo,0], (6.20) 



and 82 > —a- We omit the details here. 

This completes the proof of the proposition. □ 

Remark 6.3. We can check from formula (|6.7I) that 81 has the same sign as a±. For example, if 
a% > then 

< gQii < <5i < 50:1, 

so in particular 5i > 0. Then the density of \i\ is negative on the interval (0,<5i). Similarly, if 
o.\ < then 

|ai < <5i < |ai < 0, 

so in particular 81 < and /ii is positive on [0,/3]. 

We recall that the constants a, 81 81 and 82 in the above proposition all depend on n. From 
P^g jl -(|2^D |I . ltO ]) -([6^3 ]l and it follows that 

a = a* + 0(n- 2 / 3 ), 8 = P* + 0(«~ 2/3 ), (6.21) 
^ = L x [X-t) + L 3 t n _ m + 0(n -4/ 3)) (6 22) 

5 2 = £2(1 zl±^ »-2/3 + 0(»-*/ 3 ), (6.23) 
as n — > 00, where a* and /3* are given by (|5.36p and (I5.37p . respectively. 
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6.2 The ^-functions 

Let F\ and Fi be the Cauchy transforms of the minimizers [i\ and /Z2, i.e., 
Fj(z) = [ ^—d(J,j{x), zeC\supp(/ij)) 3 = 1 , 2 - 



2 — X 

Clearly, g'j(z) — Fj(z). In view of (|6.15[) and (|6.19l) . we have the following useful identities: 



] Vl{z)-F 1 {z) S j 



2 lw J 4t 2 (l-tfz ' 

1^ v r„\Y (z-S 2 ) 2 (z + a) 



which are valid for every z G C. 

The ^-functions are defined as follows: 



Definition 6.4. We define 



2 iw w 2i(l — f) V z 

1 X / i \ V2 



^{z) = -Vi(z)-F 2 {z) 



2 ZK ' w 2i(l — f) V 2 



(6.24) 



where the branch of the square root is taken which is positive for positive z > /3. Then £i is 
defined and analytic in C \ [0, 0\, £2 is defined and analytic in [—a, 0]. 

The ^-functions have the asymptotic behavior 

&(*) = \vl(z) - ^ + 0(z~ 2 ), fc = l,2, (6.25) 
z z 

as z — >• 00. 

6.3 The A-functions 

The functions Afc are defined as the following anti-derivatives of the ^-functions: 
Definition 6.5. We define 

Mz) = [ tk(s) da, fc = l,2, (6.26) 
Jo 

where the contour of integration does not intersect (0, 00) if k = 1 and (—00, 0) if k = 2. Then 
Ai is defined and analytic in C \ [0, 00), A2 is defined and analytic in C \ (—00, 0]. 

From (|6.24l) and (I6.26p . we have the following explicit expressions for the A-functions: 

M z)= 2 *-/-fV -^)V 2 

8t(l - 1) ^ 6>27 ^ 



pi (log - f + (js 2 - j9«) 1/2 ) - log (-#)) , 
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where the logarithm is defined with a branch cut along the positive real axis, so that for example 
log(-|) = log(f)+7ri, and 

2z + a - 4<5 2 , 2 , u/2 

X ^ = — 5771 77-( + az r 

8t(l - t) (6 28) 



Vi 



(log (z + § + (z 2 + az) 1 / 2 ) - log (§)) , 



where now the logarithm is defined with a branch cut along the negative real axis. 

Integrating (|6.25l) (or from (|6.27p - (|6.28l0 . we get the following asymptotic behavior: 

X k (z) = ^V k (z)- Pk \ogz + e k + 0(z- 1 ), fc = l,2, (6.29) 

as z — > oo, for certain constants £1,^2 that can be computed from (|6.27p - (|6.28l) . 

In what follows, we will also need the following inequalities for the A-functions. Here we write 
Xi : ±(x) for x > to denote the boundary values of Ai obtained from the upper or lower half 
plane respectively, and similarly we define A 2j ±(x) for x < 0. 



Lemma 6.6. We have 



ReAi+Or) = ReAx_(x) J °' X G [ °' ^ (6.30) 
' K ; W \>0, a; g 08, 00), 

ReA 2 . + (x)=RcA 2 ,_(x){"°' x ^^ a ^ ( 6 . 3 i) 

I > 0, x G (—00, —a), 

ImAi, + (x) = -ImAi -(x) = 717x1 ((0,x]), x > 0, (6.32) 

ImA2,+ (a;) = — ImA2.-(a;) = 7r/i2([x,0)), x < 0. (6.33) 

Proof. Using ([676]). flUJ) and (|672lj) . we find 

£1 ±( a; ) = ±7ri-(— , < x < /3, 
ax 

£2 ±(x) = ±7ri— — , —a < x < 0, 
ax 

which after integration, see (|6.26p . leads to the equalities in (|6.30[) - (|6.33p . From (|6.24[) it is also 

clear that £i(x) > for x > (3 and £2(2;) < for x < —a. This leads after integration to the 

inequalities in ([OD]) - ([675I]) . □ 

Note that as a consequence of (|6.30p and f|6 . 32[) we have 

Ai_+(x) — Ai__(x) = 27rz/xi((0, x]) = 2nipi, x > /3. 

Since np± = n\ is an integer, see (I2.15p . we then obtain 

e n(Ai 1+ (x)-Ai,_(x)) = lf ( 6 34 ) 

Similarly we have 

e n(A 2 , + (aO-A 2 ,-(aO) = ^ x < - a . (6.35) 
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7 Steepest descent analysis for Y(z) 

In this section we perform the steepest descent analysis of the RH problem 12.11 for Y, To this 
end we apply a series of explicit and invertible transformations 

Y^X^U^T^S^R 

of the RH problem. In Section [5] we will use these transformations to prove Theorems 12.71 and 

ens 



7.1 First transformation: Gaussian elimination 



In the first transformation we apply Gaussian elimination to the jump matrices of the RH problem 
for Y(z). This kind of operation was also done in |19j . Note that in QI5] there were actually two 
types of Gaussian elimination, since we needed to make a case distinction between < t < t CT n 
and t cl it < t < 1. In the present case, however, we are looking precisely at the critical time 
t = icrit, and therefore we are able to apply both kinds of Gaussian elimination simultaneously. 
We introduce a curve T r in the right-half plane passing through at angles ±y>2 with 

< ip 2 < 7T/3 

and extending to infinity and a similar curve Ti in the left-half plane with orientation as shown 
in Figure |H1 The domains enclosed by Ti around the negative real axis, and by T r around the 
positive real axis are called the global lenses [7]. 

Now we start from the solution Y of the RH problem 12.11 and define a new matrix- valued 
function X — X(z) as follows. 



Definition 7.1. Wc define 



X = Y 



X = Y 



( 1 

^1,2 



V o 

/I 








1 
1 




101,2 

1 






' 


"2,2 
li'2,1 

1 / 

0\ 



1 

«>2,1 1 I 
W2.2 / 



in the global lens 

around the positive real axis, 



in the global lens 

around the negative real axis, 



X = Y, 

where we recall the weights wj^, j, k — 1, 2 from 



elsewhere, 



(7.1) 

(7.2) 
(7.3) 



The matrix-valued function X satisfies a new RH problem, with jumps on the contour K U 
Ti U T r with the jump matrices as shown in Figure |BJ 
Thus X satisfies the following RH problem. 

RH problem 7.2. 

(1) X is analytic in C \ (M U T t U T r ). 

(2) On R U Ti U r r; we have that X + = X-Jx with jump matrices Jx as shown in Figure 



4G 




Figure 6: The curves Ti and T r , and the jump matrices Jx on R U T r U T; in the RH problem 
OforX. 



(3) j4s z — > oo, we have that 

X(z) = (I + 0(l/z))diag(z ni ,z n2 ,z-"\z-" 2 ). (7.4) 

Note that we have taken the temperature T = 1, so that by (|2.9p , (|2.10[) and (|6.ip 

ttfi,i«*,i = e-" Vl , W i, 2 h; 2 ,2 = e~ nV2 , (7.5) 
which are the non-trivial entries in the jump matrices on the real line, and 

e f (oi-oa)a! _ e -f (Vi(x)-V 2 (x)) e nKx ^ {7-6) 



w 2 .\{x) 
with a constant 



wi, 2 (s) 

e — 1=1(61-62)* _ gf (Vi(x)-V r 2 (a;))gn«:a: ) (j ^ 



1 /ai - o 2 61 — 62 1 

K = 2 — f-r)' (7 - 



which appear in the jump matrices on T r and Ti . 

The constant n depends on n, and tends to as n — > 00. Indeed we have 

k=-( ^— ^ - n -V\ (7.9) 



2 V t l-t 
on account of (|2~2"7f and (pl2"9"j) - ([230)1 . 



7.2 Second transformation: Normalization at infinity 

The next transformation is to normalize the RH problem at infinity. To this end we use the 
Afc-functions defined in (I6.26[) . 
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Figure 7: The jump matrices in the RH problem 17.41 for U . 



Definition 7.3. We define a new 4x4 matrix-valued function U = U(z) by 

U = L- n XA n , (7.10) 

where A = A(z) is given by 

A = diag (exp (Ai - |Vi) , exp (A 2 - \V 2 ) , 

exp (-Aj + \Vi) ,exp (-A 2 + \V 2 )) , (7.11) 

and 

L = dmg(e l \e l %e- l \e~ l2 ) , (7.12) 

with l x and i 2 as in 

Then U satisfies the following RH problem. 
RH problem 7.4. 

(1) U is analytic in C \ (R U Ti U r r ). 

(2) On K U T; U r r , we have that U+ = U-Jjj with jump matrices Ju as shown in Figure^ 

(3) As z — > oo, we have that 

U(z)=I + 0(\/z). (7.13) 

The asymptotic condition of U in (|7.13[) follows from (|6.29[) and the jump matrices Jjj in 
Figure [7] follow from (|7.5[) ~ (|7.8[) and straightforward calculations. 
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By (|6.30[) - (I6.35[) the jump matrix Jjj takes the following form on (—a, f3), 



Ju 



Ju 





2nA 2 ,+ 



/I 

e 

\0 

/ a 2nXi, + q 

1 







1 



= 2raA 2 



1 



o 



\ 

■-/ 
°\ 




V 



on (—a, 0) 



on (0,0) 



(7.14) 



(7.15) 



and the following form on the rest of the real line 

/ + e -2„RcA 1 , + £i3 



Ju 
Ju 



on (— oo, —a), 
on (j8,oo), 



which tend to the identity matrix as n — > oo at an exponential rate because of the strict inequal- 
ities in (I6.30[) and (|6.31[) . Recall that E^j denotes the elementary matrix with 1 at position 
and zero elsewhere. 

7.3 Third transformation: Opening of local lenses 

In the next transformation U n- T we open a local lens around each of the intervals [—a, 0] and 
[0,/?]. 

The local lenses around [0,(3] and [— a, 0] are illustrated in Figure |S] The lips of the lens 
around [0,(3] are denoted by Ei, and the lips of the lens around [— a, 0] are denoted by £2- We 
make sure that Si and S2 do not intersect with r r and T;, except at the origin. 

The transformation U H> T is based on the standard factorization (we use + Aj ; _ = 0) 







1 

2nA, 



1 

,2nA^ 



1 

-1 



1 

3 2nAj,+ 



that can be applied to the non- trivial 2x2 blocks in (17.141) and (|7. 15|) . 
Definition 7.5. We define 



T(z) = 



T{z) 



U(z)(I- 
U(z)(I- 

'U(z)(l- 



2nAi 



£3,1) in upper lens region around [0, 0\, 
£3,1) in lower lens region around [O,0\, 



e £'4,2) in upper lens region around [— a, 0], 



U (z)(I + e 2 ™ 2 £4,2) in lower lens region around [—a, 0]. 
T(z) — U(z), outside the lenses. 

Then the matrix- valued function T is defined and analytic in C \ Sy, where 

s T = mu r ; u r r u Si u s 2 

and it satisfies the following RH problem. 
RH problem 7.6. 

(1) T is analytic in C \ S^. 



(7.16) 

(7.17) 
(7.18) 
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Figure 8: The jump contour Et in the RH problem 17.61 for T. 



(2) T satisfies the jump relation T + = T-Jt on E^ with jump matrices: 



Jt 



Jt 



Jt 



Jt 




on E; 



on E 



on (—a, 0), 



on (0,(3) 



■2: 



1- 



and 



Jt — Jui 



on Ti U r r U (—oo, —a) U (/3, oo). 



(3) As z — > oo, we have that 



T{z) = I + 0{l/z). 



See Figure [5] for the jump contour Et and Figure |H] for all jumps Jt in a neighborhood of 
the origin. 

7.4 Large n behavior 

Now we take a closer look at the jump matrices Jt- The jump matrix is constant on [—a, 0] and 
on [0,(3]. On all parts of E^ \ [— a, 0] the jump matrix is equal to the identity matrix plus one 
or two non-zero off-diagonal entries. Ideally, we would like to have that all off-diagonal entries 
tend to as n — > oo. However, we cannot hope for this in a neighborhood of 0, due to the 
fact that we modified the equilibrium problem near 0. The exceptional neighborhood of is 
shrinking as n increases at a rate of 0(n~ 2 / 3 ). Later we will construct a local parametrix in a 
larger neighborhood of of radius O^™ 1 / 3 ). 
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Figure 9: The jump matrices Jt in the RH problem l7.6l for T in a neighborhood of z — 0. A local 
parametrix in this neighborhood is built by means of the RH problem 12.21 for M as in Figure [5] 



We show here that outside this shrinking neighborhood the jump matrices Jt on \ [—a, 0\ 
do indeed tend to the identity matrix as n — > oo. 

The functions Ai and A2 that appear in Jt depend on n. They have limits as n — > 00, which 
we denote by AJ and A|. Since f3 /3* , a ->• a* , 5i,a -> (see l|6.2ip - (|6.23p ) and pi j2 -> p* 2 as 
n -> 00, we have by (|6.24l) and (I6.26j) , 

= / fc = l,2, (7.19) 

Jo 

where 

where the branches of the square roots are taken with are positive for large positive z. 
We also have 



- p{ (log (* - £ + <2 2 - ,T=) 1/J ) - log (-£)), (7-21) 

- (log (z + f + (z 2 + a-z)" 2 ) - log (if )) , (7.22) 
with the appropriate branches of the logarithms. 
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Proposition 7.7. There exists a constant C > such that for every large n, we have 

\\ k {z) - \* k (z)\ <Cn- 2 / 3 max(\z\,\z\^ 2 ), z £ C, (7.23) 

for k = 1,2. 

Proof. We have by ([B~2"H) . 

<5i(z - /3) 1 / 2 zV2(( z _ 0)1/2 _ ( z _ /3*)l/2) 



2t(l-t)z 1 /2 2t(l-t) 

s^z-py/ 2 {p-p^z 1 / 2 



2t{l - t)zV 2 2t(l - t){(z - /3)!/2 + [z - /3*)i/2) 

h{z-fl) 1/2 P~P* z^ 2 f (z-/3)V 2 \ " 
2t(l-t)z 1 / 2 2t(l - t) (z - p*) 1 / 2 \ (z-P*) 1 / 2 ) 



(7.24) 



Note that f _a,\ 1/2 has non-negative real part, so that 



(Z-/3*) 1 / 2 

l 

< 1 



1+ (^-/3) 1/2 



for every z. Then since = 0{n~ 2 / z ) and P — P* = 0(n~ 2 / 3 ), we can estimate (17.24)) as 

- (*) = OCn- 2 / 3 ^!" 1 / 2 ) + 0(n- 2 / 3 |z - rr 1 / 2 ) + 0(n- 2 / 3 ), as n -> oo, 

uniformly for z € C. After integration from to z, see (|6.26p . we find (|7.23p with k — 1. The 
proof for fc = 2 is similar. □ 

7.5 Estimate on local lenses 

We need an estimate for ReAi on Ei outside of the shrinking disk of radius n -1 / 3 , 

L>(0,n- 1/3 ) = {z e C | |z| < n- 1 / 3 } 

around and a fixed disk D(P*, e) around P* . The following lemma gives such an estimate, 
together with a similar estimate for ReA2 on £2. 

Lemma 7.8. (a) Let e > 0. Then there is a constant c\ > such that for every large enough 
n, we have 

ReAi(z) < -c in - 1/2 , ze Si \(D(0,n- 1/3 )U D{P*,e)). (7.25) 

(b) Let e > 0. Then there is a constant C2 > smc/i £/ia£ for every large enough n, we have 

ReA 2 (z) < -c 2 rr 1/2 , z6S 2 \(J)(0,^ 1/3 )UD(-« , ,f)). (7.26) 

Proof. For A| it is easy to show from (|7.19[) and (|7.20j) that Re A|(z) < in the strip < Rez < 
P*, Imz ^ 0, and 

ReAJ(z) < -c' im in(|z| 3 / 2 ,|z-/3*| 3 / 2 ), ze^ 

for some constant c[ > 0. Then by (|7.23[) we have for another constant c" > and for n large 
enough. 

ReAi(z) = ReAJ(z) + Re(Ai(z) - A^(z)) 

< -c'[\zfl 2 + C'n- 2 '^' 2 : z e Ei \ D(P*,e) 

with C" > 0. This leads to (|7.25l) with a suitable constant c% > 0. 

The proof of part (b) is similar. □ 



52 



7.6 Estimate on global lenses 

We start by being more precise on the location of the global lenses T; and T r . 

Lemma 7.9. We can (and do) choose the global lenses T r in the right-half plane and Ti in the 
left-half plane such that for some constants c r ,ci > 0, 

Re(AJ(z)- X* 2 (z)) < -c r min(|z|,|z| 3/2 ), forz£T r , (7.27) 

and 

Re(AJ(z) - X* 2 (z)) > cimm(\z\,\z\ 3 /% for z e T,. (7.28) 
Proof. We consider the curve 

C: Re(A*(z) - X 2 (z)) = 0. 
By definition A*(0) = X 2 (0) = and so the curve C contains 0. By (I7.20[) we have as z — > 0, 

(*) = i0r/ /2 + °(- 3/2 ) = ^rr^ ^ 1/2 + 0(z3/2) < ±Imz > °> 

^ 2[Z) - 2t(l-t) Z +U[Z >~ i3/4 (1 _ i)3 /4 Z + J> 

where we used (|5.36|) and (|5.37l) . After integration, we find 



(7.29) 



and so 



A 1 (z)-A 2 (z) = 3^/4^ _ ^3/4 ~ z +0{z 0/ ), ±lmz>0, 

as z -> 0. This implies that the curve C makes angles 

± - + - arctan ' 



3 3 V(P2) 1/4 

with the positive real axis. If p* = p 2 = 1/2, then this is a right angle, and then in fact C 
coincides with the imaginary axis. 

As z -> oo, we have by (|7T2Tj) and (jTT^j) . 

(a* + 0*) 

Xl(z) - X* 2 (z) = - \ t ^_ t ^ + OQog \z\), 

which implies that C gets more parallel to the imaginary axis as z — > oo in the sense that 
arg z — >■ ± § as z — > oo with zeC. 

Thus C divides the plane into a part on the right where Re (X*(z) — X 2 (z)) < and a part 
on the left where Re (X*(z) — X 2 (z)) > 0. We can then take a smooth curve T r in the region to 
the right of C, making angles ±</?2 with < tp 2 < 7r/3 at the origin, and extending to infinity in 
such a way that (|7.2T|) holds. 

Similarly we can take T; to the left of C such that (|7.28[) holds. □ 
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Now we can make the necessary estimates for Re (Ai — A 2 ± n n z) on r r and Ti outside of 
the shrinking disk Z?(0, n" 1 / 3 ) of radius n -1 / 3 around 0. Here we write k = K n to emphasize its 
dependence on n, see (|7.8|l and (|7.9D . 

Lemma 7.10. (a) There is a constant c\ > swc/i that for every large enough n, we have 

Re (Ai(z)- A 2 (z)±kz) < -ci min(|z|, |z| 3/2 ) < -cin~ 1/2 , z G r r \L>(0, n _1/3 ). (7.30) 

(b) There is a constant c 2 > such that for every large enough n, we have 

Re (Ai(z) - A 2 (z) ± kz) > c 2 min(|z|, |z| 3/2 ) > c 2 n~ 1/2 , z G T ; \ D(0, n~ 1/3 ). (7.31) 

Proof. This follows by combining (|7.23[l with the estimates (|7.27p - (|7.28[) and the fact that « = 
K n = 0(n~ 2 / 3 ), see (I7.9p . We will not give details as the proof is similar to the proof of Lemma 

rrsi □ 

As a result of Lemmas 17.81 and 17.101 we have that outside the shrinking disk D(0,n -1 / 3 ) the 
jump matrices Jt on the local and global lenses tend to the identity matrix. More precisely: 

Corollary 7.11. (a) For every e > there is a constant c\ > such that 

J T (z) = I + C>(e- Cinl/2 ) as n->oo, 



uniformly for z £ (Si U E 2 U T ; U T r ) \ (D(0, n" 1 / 3 ) U D(-a* , e) U D(f3, e)) . 

(b) There is a constant c 2 > suc/i that 

J T (z) = I + O (e- C2 "l z l) asn^oo, 

uniformly for z £ T; U T r luii/i |z| > 1. 

(c) For every e > i/iere is a constant C3 > such that 

J T {z) = I + (e- C3 "l 2 l 2 ) asn^ao, 



uniformly for z 6 (—00, —a* — s] U [/J* + £, 00). 

Proof. Part (a) follows from the estimates in Lemmas l7.8l and r7.10l which apply to the off-diagonal 
elements in the jump matrix Jt on the lips of the local and global lenses. 
Part (b) also follows directly from Lemma 17. 101 

Part (c) follows from the variational inequalities in (|6.30[) and (|6.31l) that in fact can be 
strengthened to 

Re Ai i+ (x) = ReAi,_(ir) > c^x 2 , x>f3 + e, 
Pe A 2j +(x) = ReA 2 ,_(x) > c^x 2 , x < —a — e, 

for some constant C3 > independent of n. □ 
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7.7 Global parametrix 

In this section we construct a global parametrix P(°°>(z) for T(z). The matrix- valued function 
P(°°'(z) satisfies the following RH problem, which we obtain from the RH problem 17.61 for T by 
ignoring the jumps on the local and global lenses and on (—00, —a) and (f3, 00). That is, we keep 
the jump on [—a, 0\ only. 

RH problem 7.12. 

(1) P(°°)(z) is analytic for z e C\([-a, /?]). 

(2) p(°°) satisfies the jumps 



(00) , 
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/or x € (—a, 0). 



(3) As z — ¥ 00, we have that 



P ico \z) = I + 0(l/z). 



A solution of the RH problem 17.121 is 



P^'tz) = - 
y ' 2 



where 



(1 i 0\ 
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(7.32) 



(7.33) 



and where we choose the branches of the 1/4 powers that are positive and real for large enough 
real z. 



7.8 Local parametrices around the non-critical endpoints 

On a small but fixed neighborhood around each of the non-critical endpoints —a and /?, we 
construct local parametrices p(~ Q ) and P^ out of Airy functions. We can in fact take n- 
indcpendent disks D{fi*,E\) and D(— a*, £2) for certain £i,£2 > such that P^ and 
have the same jumps as T has in the respective disks and such that 

P {l3 \z) = (I + 0{n- 1 ))P {aa) (z), uniformly for \z-fi*\= £ X) 
P ( ~ a \z) = {I + O(n- 1 ))p(° o) (z), uniformly for \z + a* \ = e 2 , ' '"' J) 

as n — > 00. 

The construction with Airy functions is well-known in the literature and we do not give 
details, cf. also Section |3~71 of this paper. 
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7.9 Local parametrix around the origin 

7.9.1 Statement of the local RH problem 

In this section we construct a local parametrix around the origin, with the help of the model RH 
problem for M(£). We want to solve the following RH problem. 

RH problem 7.13. We look for P^ satisfying the following: 

(1) P(°(» is analytic for z € D(0 : n 1 / 3 ) \ Tjt, where D(0,n 1 / 3 ) denotes the disk of radius 
n" 1 / 3 around 0. 

(2) P(°) satisfies the jumps 

p|0) = p(0) m ^ n i J( 0jn -l/3) ) 

where Jt is the jump matrix in the RH problem for T. 

(3) As n — > oo, we have that 

P (0) (z) = (l + N{z) + 0(n- 1/3 )^j P (oo) (z) uniformly for |z| = nT 1 /3 (7.35) 

where 

N(z) — 0(1) as n — > oo uniformly for \z\ — n -1 / 3 . (7.36) 

The matching condition in (3) is posed on the circle \z\ = n^ 1 / 3 which is shrinking as n 
increases. It is different from the usual matching condition which requires 

P(°)(z) = (I + 0(n- s )) P ico) (z) asTWoo (7.37) 

for some S > 0. It turns out that in general we cannot achieve (I7.37p . Only in case n = we 
can achieve (|7.37l) on a circle of fixed radius. For k ^ 0, we have to control the terms e ±nKZ that 
appear in the jump matrix Jt on Ti and T r . These terms remain bounded if \z\ = 0{n-^ 3 ) 
which (partly) explains the radius n" 1 ! 3 of the shrinking disk. 

An essential issue for the further analysis is that the 4x4 matrix N(z) from (|7.35[) - (|7.36[) is 
nilpotent of degree two, i.e., N 2 (z) — 0, and even more 

N(zx)N(z2) = for any z±, z-i different from 0. 

The matrix-valued function N(z) is also analytic in a punctured neighborhood of with a simple 
pole at 0, see the explicit formula in Section f7.9.4l See also [2TJ [38] for a similar feature in the 
RH analysis. 

7.9.2 Basic idea for the construction of the parametrix 

We will construct P^ ' in the form 

P<°)(z) = P (0) (z) diag (e nXl< ^ z \ e "(^ e -«Ai W ) e -n(A a (*)-««)^ _ 

Then in order that P^ has the jump matrices Jt, the matrix-valued function p(°) should have 
constant jumps on each part of St H -D(0, n" 1 / 3 ) and these constant jumps are exactly the same 
as the jumps in the RH problem l2.2l for M, cf. the jump matrices in Figures [5] and [SJ The model 
RH problem 12.21 depends on parameters r%, f2, si, S2- It follows that for any choice of these 
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parameters, any conformal map f n {z) that maps the contours £t (1D(0, n - 1 / 3 ) into the ten rays 
Uj=o ^-j anc ^ an y analytic prefactor E n (z), the definition 

Pt°\z) = E n (z)M {f n (z);r 1 ,r 2 ,s 1 ,s 2 ) 

x diag (e nXll - z \ e n ( x ^( z )+ KZ ) ^ e -» lA i(«) ; e -n(A 2 («)-«*) 



will give us P^ that satisfies the jump conditions in the RH problem l7.13l for PW . The conformal 
map /„ and the analytic prefactor E n (z) should then be chosen in order to have the asymptotic 
condition (|7.35|) . 

However, there is not enough freedom to achieve this. Since the jumps in the model RH 
problem for M do not depend on the parameters r\,r%, s\, s 2 , we also let these depend on z and 
on n, and we put 

P(°>0) = E n {z)M (f n ( z ); n,„(z), ra, n («), ai,»(«), «a,„(*)) 

diag (e nAl ^, e Ii ( A 2( 2 )+ K2 ) i e -n^iW i e -«( A 2(z)-K2)^ ^ 

where n n (2), r 2,n(z), si : „(z), s 2 ,„(z) depend analytically on z. Now we have these functions at 
our disposal, as well as E n (z) and f n (z), to achieve the condition (3) in the RH problem 17. 131 for 
P(°). 

It turns out that the right choice for these functions takes the form 

f n (z) = n 2/3 f(z), r j>n (z) = rj(z), s J> (z) = n 2/3 Sj(z), j = 1,2, 

for certain analytic f{z), J*i(z), r 2 (z) that do not depend on n, and analytic s\(z), s 2 (z) that 
still depend on n but only in a mild way. We next describe these functions. 

7.9.3 Auxiliary functions 

Recall the functions Ai and A2 from Section 16.31 and their limits A^ and X 2 as n — > 00 from 
Section O 

We define the following auxiliary functions. 

Definition 7.14. (a) We define 

2/3 



Hz) = -( Pl r 1/6 (§*!(*))' 



(7.38) 



which is a conformal map in a neighborhood of the origin, see also (|7.43[) . 

(b) We define functions 

r 1 =r 1 (z) = (pl) 1 '\ r 2 (z) = iifc^Mf) ±I m2 > , (7.39) 

A 1 (Z) 

where 7*1 is a constant and r 2 (z) is an analytic function near 0, independent from n. 

(c) Finally, we define the n-dependent functions 

1 Xi(z) - Xt(z) 1 X 2 (z) - X* 2 (z) 

Sl(2)= 2 Ff{zW 2 ' 2i ) = 2 /(z) 1 ^ ' (7 ' 40) 
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Note that the functions are defined such that 



A 2 (z) = ^ 2 (z)(/(z)) 3 / 2 + 2 S2 (z)(/(z)) 1 /2. 



(7.41) 
(7.42) 



Lemma 7.15. There is r > such that the junctions f , r 2; Si and s 2 are analytic in the disk 
D(0,r) — {z | \z\ < r}. In addition the following hold. 



(a) The function f(z) is real for real z S D(0,r) ; and 



f{z) 



+ 0(z 2 ) 



as z — > 0. 



(b) The function r 2 (z) is real and positive for real z 6 D(0,r) and 

ra(0) = ij>l) 1/A - 



(7.43) 



(7.44) 



(c) The functions s\ and s 2 depend on n in such a way that n 2 ^s\ and n 2 / 3 s 2 have limits as 
n — > oo. The limiting functions are analytic in D(0,r) and satisfy 



lim n 2/3 si(0) 

n— >oo 

lim n 2/3 s 2 (0) 



(Pi 



,1/4 



2(\M + Vp!) 



(^) 1/4 



(7.45) 



2(^pT + ypI) 



where L$ and Lq are given in (|2.31|) — f|2 .32[) . 

Proof. Parts (a) and (b) follow from the definitions of AJ and A 2 in (|7TTg|) - (|7T2lI)) . see also (|7T2^) . 
For part (c) we note that by (j7^H|) . (j6~Ml) - (joT2l)j) . and (fTTT^l) (f7T2l7|) . 



lim n 2/3 si(z) 



lim 



™ 2/3 / (tii(s)-&(s))ds 



2(-/(z))V2 „^oo 

hm™ [n 2 / 3 (/?-/?*)] / • ^/ 2 



2t(l-t)(-/(z))V2 



s-/3* 



lim„_ 



2/3 <5il 



4t(l-t)(-/(«))Va J 



s-/3* 



1/2 



(is, 



(7.46) 



which by (|7.43p is indeed an analytic function in a neighborhood of 0. Recall that j3 and Si 
depend on n, and that the limits of n 2 / 3 (/3 — /?*) and n 2 / 3 o"i as n — > oo exist. 

To evaluate (|7.46p for z = we only need to consider the second term on the right-hand side 
of ([LHP . By (jB~m (fP31 and (jT^Sl) we then obtain 

2 /3 / n x Li(l-t)+L 3 * /Ta- 
llin n / J si(0) = ^ qM/1 irTTFVP 



4t 3 / 4 (l - i) 3 / 4 

which can be rewritten to the formula given in (|7.45[) using (15 -3T|) . (|2.26l) . (12.27[) and (|2.31|) . 
The statements in part (c) dealing with s 2 follow in a similar way. □ 

Observe that n in ([73^1) . r 2 (0) in (TTi4]) and the limits of n 2 / 3 si(0) and n 2 / 3 s 2 (0) from (f7T^5|) 
correspond precisely to the values in (|2.50[) (|2.51l) . 
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7.9.4 Definition of parametrix 

We finally define the local parametrix P'°) near the origin as follows. 
Definition 7.16. We define 



pW(z) = £;„(z)A/(n 2 / 3 /(z);r 1 (z),r 2 (z),n 2 / 3 Sl (z),n 2 / 3 S2 (z) 



x dia; 



iag (e" Al 



(z) n(\ 2 (z)+K,z) -nAi(z) —n(X 2 (z)-Kz) 



(7.47) 



where n is given in (|7.8p and 

E n (z) = p(°°\z)di ag (l,e- 



«,l,e- 



/l £ 0^ 

1 1 -£ 
} V2 » 1 
\0 -£ 1 J 

x diag ((-n 2 / 3 /(*)) 1/4 , (n 2/3 /(^)) 1/4 , {-n^f{z))-y\ (n^f(z))-^) . (7.48) 

Remark 7.17. Note that the parameters ^(z), rj 2 / 3 si(z), n 2 / 3 S2(z) of M in (|7.47p are real for 
real z, but in general non-real if z is not real. We proved in Theorem l2 . 31 that the RH problem for 
M is solvable, and thus that M(£; r%, r<z, Si, S2) exists for real parameters (and r\ > 0, r 2 > 0). 
By perturbation arguments it can then be shown that the RH problem for M is also solvable for 
parameters that are sufficiently close to the real line. 

We are interested in z G Z^O,™™ 1 / 3 ) and for such z, the values of r 2 (z), n 2 / 3 s\(z), n 2 / 3 S2(z) 
come arbitrarily close to the real axis as n — >■ oo. Therefore for large enough n, M (£; ri(z), T2{z) , 
n 2 / 3 si(z), n 2 / 3 S2(z) s ) exists for every z G -D(0, n -1 / 3 ) and the local parametrix (|7.47p is well- 
defined. 

Also the asymptotic condition in the RH problem 12.21 for M will be valid uniformly for 

z G D^.n^l 3 ). 

Lemma 7.18. The prefactor E n (z) in (|7.48p is analytic in a neighborhood of the origin. 
Proof. By using the expression for p(°°>(z) in (I7.32p . we can rewrite (17.481) as 



E n (z) = diag(l,e~ 



! ,l,e~ 



>V2 



I 1 





% 


°\ 





i 





—i 


i 





1 





v° 


— i 





1/ 



x diag ( 7l (z)(-n 2 / 3 /(z)) 1 /4 j7 -i( z )(„2/3 /(z)) i/4 j 

7 r 1 (z)(-n 2 / 3 /(z))- 1 /4 ;72(z)(n 2/3 /(z)r i/4^ (749) 

with 71 and 72 as in (|7.33p . Here we used that the matrix diag(l, e~ nKZ , 1, e~ nKZ ) commutes 
with all the other matrices in (T71g)) . From ([TI9")) . (1713)) and ([7T55]l . we then see that P„(z) is 
indeed analytic in a neighborhood of the origin. □ 

As discussed in Section 17.9.21 the matrix- valued function PW defined in (|7.47p satisfies the 
jump condition in the RH problem 17.131 [We modify the contours St if necessary, in such a 
way that / maps St H D(0, n -1 / 3 ) into U^=o ^V] ^ then remains to show that p(°) satisfies the 
matching condition (3) in the RH problem 17.131 
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Lemma 7.19. P'°) defined in (17.47[) satisfies the matching condition (3) in the RH problem 

[7731 

Proof. The asymptotics (|2.37p for M(£) can be rewritten as 



m(o = dia g ((-c)- 1/4 , r 1/4 , (-o 1/4 , c 1/4 )^ 



/l 





—i 


0> \ 





1 





j 


—i 





1 





\o 


i 





V 



92(C) je 9i(C) ;e e 2 (C) 



(7.50) 



where 



Mi = - 



(1 i 0\ 

1 -i 

1 1 
\0 -i 1 / 



/e^ 7 "/ 4 0\ 

10 



\ 0/ 



Ml 



/o 








°\ 


(1 





—i 



















1 





/' 













—i 





1 





\o 








V 


\o 


i 





V 



as ±Im C > 0, and M 2 takes the following form: 



Mo 



1 -i 

1 1 
\0 -i 1 



\ 


(* 


* 







/I 





—i 


°\ 




* 


* 











1 















* 


* 


— i 





1 





/ 







* 




\o 


i 





1/ 



with * denoting certain unimportant entries. 

Recall that M(£) as well as Mi depend on ri,r2,s%,S2- Also #j(C) depends on Tj and Sj, for 
j = 1,2, see (|2.38l) - (|2.39p . Using the notation 9j((;rj, sj) to denote the dependence on fj and 
Sj, we have 

^■(n 2 / 3 /W;r-,(z),n 2 / 3 S ,(z)) = nA,(z) 

because of ([2~3g[t - ([2~39T> and ([7H ]) -([732" ]1 . 

Then we have for z on the circle \z\ = n~ 1/3 , by ([7I7[) . ([715)1 and ([730"[i . 



(7.51) 



p(°)(^)(p(«))- 1 ( 5: ) 

= P(°°)(z)diag(l, 



\l,e-" KZ ) 7 + 



Mi (2) 



n!/3/i/2( z ) 



x diag(l, e" KZ , 1, e™ KZ )(p( °))- 1 (z), (7.52) 



as the exponential factors involving 9j and Aj cancel due to (|7.51[) . In (|7.52[) we have written 
Mi(z) for 

M 1 (r 1 (z),r 2 (z),n 2 ^ Sl (z),n 2 ^ S2 {z)) 

by a slight abuse of notation. 

Using the formula ([7521 for p(°°\z) in ([7321 . we obtain 



pW(z)(p( °))- 1 (z) = 7 + 2V(s) + O^' 1 / 3 ), for |z| = yi" 1 / 3 , 



(7.53) 



GO 



with 



N(z) 



(1 i 0\ 

1 -i 

1 1 
\0 -i 1 / 



diag(e =F7ri / 4 7i(z) e -" KZ 7 2 " 1 (z) 0) 



2n l/3/l/2( z ) 

xMi(z)diag(0 e^^O) e" KZ 7 2 ~ 1 (^)) 



/l 





—i 


°\ 





1 





/ 


— i 





1 





\o 







V 



(7.54) 



where 



M^z) = M 1 (r 1 (0),r2(2) I n 2 / 3 s 1 (0),n 2 / 3 s 2 ^)) 



is analytic. 

Then by the definition of M 1 and 4731) . 



iV(z) 



(l i 0\ 
1 1 -i 

2n 1 /3/i/2(^) j 1 
\0 -i 1 / 

/0 ±7 2 (z)c(z) ie =F W4 e n^ 7l ( z ) 7 -i( 2 ) d ( z )^ 

ie^ 4 e- nKZ -f 1 (z)j 2 ; 1 (z)d(z) i^ 2 (z)c(z) 



\0 / 



(I Q -i Q\ 
1 i 
-i 1 
\ i 1/ 

(7.55) 



with d(z), c(z), c(z) given by (|2.41[) - (|2.43p with n, r 2 , Si, S2 replaced by ri(z), r 2 (z), n, 2 / 3 si(z), 
n 2 / 3 S2(2;), respectively. These functions are analytic and remain bounded as n — > oo. 
The combinations 



71 (z) 



1 



72(Z)/V2( Z )' 7 |( Z )/l/2 (z) ' 



=Flmz > 0, 



are analytic in a punctured neighborhood of with a simple pole at 0. On the circle with radius 
n ' 3 they all grow like 0{n^ 3 ) as n — > oo. Since in (|7.55[) there is also a factor n 1 / 3 we find 
that indeed 

N(z) = 0(1), as n — ► oo uniformly for |z| = n~ 1//3 , 



and this proves the lemma. 



□ 



From the proof of the lemma, we also find that N(z) is analytic with a simple pole at 0. From 
([734]) or ([735)1 it is also clear that 



N{z\)N{z2) — 0, for any zi, z^ different from 0. 
7.10 Fourth transformation 



(7.56) 



Using the global parametrix p(°°) and the local parametrices a \ P^ and P(°\ we define 
the fourth transformation T i— > S as follows. 
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Definition 7.20. We define 



{r(p(0))-i 5 insi de the disk Z^O,™ -1 / 3 ) around 0, 

T(P(-"))- 1 , in a fixed disk around - a*, 

T(PW)- 1 , in a fixed disk around /?*, ( ' ^ 

r(pM)-i ) elsewhere. 

Then S 1 is defined and analytic outside of £t and the three disks around 0, —a* and /?*, with 
an analytic continuation across those parts of St where the jumps of the parametrices coincide 
with those of T. What remains are the jumps on a contour E5 that consists of the three circles 
around 0, —a*, and /3, the parts of the real intervals (—00, —a*] and [/?*, 00) outside of the disks, 
and the lips of the local and global lenses outside of the disks. 

The circles are oriented clockwise. Then the RH problem for S is 

RH problem 7.21. 

(1) S is analytic in C \ E5. 

(2) S satisfies the jump relation S+ = S-Js on E5 with jump matrices: 



{p(0)(pM)-l on thf, boundary of D(0, n" 1 / 3 ) 

p(- Q )(p( oc ))- 1 on ^/j e boundary of the disk around — a*, 
ir\i ( u_i (7.58) 

pyp) (p(°°)) t on i/j e boundary of the disk around p*, 

p(-°°)j T (pM)-i elsewhere on E s . 



(3) yls z — > oo 7 we /iaue i/iai 

5(2) = I + 0(l/z). 

The jump matrix J5 is not close to the identity matrix on the circle around 0, since by (|7.58[) 
and (17.351) we have 

J s = I + N(z) + 0(n- 1/3 ) asrwoo (7.59) 

uniformly for \z\ — n -1 / 3 , with N(z) — 0(1). 

The other jump matrices, however, are close to the identity matrix as n gets large. Indeed 
from (|7.58[) and the matching condition (|7.34[) we find that 

J s = I + 0(n- 1 ) asTWoo (7.60) 

for z on the fixed circles around —a* and j3* . 

From Corollary 17.111 and (I7.58P , it follows that Js has similar estimates as Jt on the rest of 
E5. Namely, for certain ci, C2, C3 > we have 

J s (z)^I + 0(e- c ^ 1/2 ) (7.61) 
on the lips of the local and global lenses outside the disks, 

Js(z) = 7 + 0(e- C2n|z| ) (7.62) 
on the lips of the global lenses with \z\ > 1, and 

J s {z) = I + 0{e- c ^) (7.63) 
on the real line outside of the disks around —a* and j3*. 
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7.11 Final transformation 

The jump matrix (|7.59[) on the shrinking circle \z\ = rt -1 / 3 around the origin does not tend to 
the identity matrix as n — > oo. The final transformation S i— > R serves to resolve this issue. An 
important role will be played by the special structure of the matrix N in (|7.54|) , see also [2TJ [35] . 

Recall that N(z) is analytic for z in a punctured neighborhood of with a simple pole at 
z = 0. Define 

N = Res 2=0 N(z) (7.64) 
as the residue matrix. Then we have the splitting 

z 



N(z) = N(z) - 



where N(z) — Nq/z is analytic inside the disk D(0 : n - 1 / 3 ) and Nq/z is analytic outside. Note 
also that because of (|7.56p we have 

N(z)N o = 0, iV 2 (z) = 0, Nq=0. (7.65) 

The transformation S i— > R is now defined as follows. 

Definition 7.22. We define 

S(z) (i + N{z) - — V for 2 £ D(0, n" 1 / 3 ) \ E 5 , 

I / - — J , for z e C \ (13(0, n-V3) y ^ 

where and ATq are given in (|7.54l) and (17.64[) . respectively. 

Then R is defined and analytic in C \ where = S5 and R satisfies a RH problem of 
the following form. 

RH problem 7.23. 

(1) R is analytic in C \ 

(2) i? satisfies the jumps R + = R-Jr on 

(3) R(z) = I + 0(1/ z) as z -)• 00. 

The jump matrix Jr for |z| = n -1 / 3 is by (|7.66[) (recall that we use the clockwise orientation) 

J R {z) = (l + N{z)-^y J s (z) (/-^ 

= (/ - N(z) + ^) (/ + iV (z) + O^- 1 / 3 )) (/ - ^) 

where the calculation of the inverse oi I + N(z) — Nq/ z was done with the help of (|7.65D . Then 
if we expand this product of three matrices, many terms cancel due to (|7.65[) . What remains is 
simply 

J R (z)=I + 0(n- 1 ^), \z\=n- 1 '\ (7.67) 



where we also use the fact that N(z) as well as Nq/z are 0(1) for \z 



= n-V3. 
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The transformation (|7.66p does not change the jump matrices on the other parts of in 
an essential way. Thus Jr tends to the identity matrix on these parts as well, with a rate of 
convergence that is the same as that for J$, see (I7.60l) - (|7.63p . 

We have now achieved the goal of the steepest descent analysis of the RH problem. The jump 
matrices for R tend to the identity matrix as n — > oo, uniformly on S^j as well as in L 2 (Yir). By 
standard arguments, see |16j . (and also [5] for the case of contours that are varying with n), we 
conclude that 

R ^ = I + {n^m) (? - 68) 
as ti — y oo , uniformly for z £ C \ E^. 

The estimate (|7.68p is the main outcome of our Deift-Zhou steepest descent analysis for the 
RH problem 12. ll for Y. We will use it to prove the main theorems in the next section. 



8 Proofs of Theorems [2T] and [2341 
8.1 Proof of Theorem I27TI 

In this section we prove Theorem 12.71 by following the subsequent transformations of the RH 
problem. We start with the expression (I2.14j) for the correlation kernel K n (x,y): 



K n {x,y) 



2ni(x — y) 



(0 w 2tl (y) w 2 ,2{y))Y+ 1 (y)Y+{x) 



wi i2 (x) 


V o J 



(8.1) 



We first assume that both x and y are positive and close to 0. From the first transformation 
F^lin ([7TT |l -(ff3 |l . it follows that 



1 



K n (x,y) = — -(0 w 2A (y) 0) X+\y)X + (x) 



2ixi(x — y) 

Applying the second transformation X i-> U (|7.10p to (18.21) , we have 

e ™(4T0T)(2/ 2 -* 2 )+ fc fe-*)) 



/ 101,1 (x)\ 



V o J 



(8.2) 



K n (x,y) = 



2m(x — y) 



(0 e-™ Al -+^ 0)U+ 1 (y)U+(x) 



/ e -nAi, + (a)\ 




V o ) 



where fc = + 2( ^ t) . This constant depends on n and from (j2~29l) . (1230)) and (12^27) it follows 



that 



fc* := hm k = -21 i 6 i 



2t 2(1 -t) 2t 2(1 -t) 



(8.3) 



with the convergence rate 0{n 2 / 3 ). 



G4 



By the third transformation U > T in (|7.16p - (|7.18[) and our assumption on x, y, it is readily 
seen that 



K n (x,y) = - — — 3 (_ e nAi,+ (3/) o g-nAi.+ fe) 



2iri(x — y) 



x T+\y)T+(x) 



/ e -nA li + (i)\ 



; nAi, + (» 

V o J 



(8.4) 



For z e Z?(0, n~ 1/3 ), we have from (17371) and (jTM)) that 

T(z) = S(z)P ( - a \z) = R{z) \I + N(z) - 
= R(z) - N(z) + P(°)(z). 
Substituting ([831) and (j7T47j) into (f8T4|) . it then follows that 



iVo 

a; 



Pf)(«) 



(8.5) 



K„(x,y) 



2iri(x — y) 

x(-l 1 0)M+ 1 (« 2/3 /(2/);ri(y),r 2 (y),n 2 / 3 Sl (y),n 2 / 3 S2 ( 2/ )) 



x E^ 1 (y) ( J + iV(y) - ^° ) RT l {y)R(x) ( 7 - N(x) + ^ ) E n (x 



x M + (n 2 / 3 /(x); ri(a;), r 2 (x), n 2 / 3 si(x), n 2 / 3 s 2 (a;)) 




1 

w 



(8.6) 



Now we fix u, v > and take 



l 2 /3' 



cn 



2/3 ' 



(8.7) 



where c = , 1 with t in (|2.27D . Then for n large enough, a; and y are inside the disk 

y/t(l— t) 

D(0, rt -1 / 3 ). It is then readily seen that 



lim e'^ 

n— >oo 



and by (J733J), 



n 2/3 f(x) -> it, n 
as n — > oo. From (18.71) we also find that 



"(i^yfe 2 - 2 )) = 

/3 /(2/) -> «, 



(8.8) 



r 2 (a;)->r 2 , r 2 (y)->r 2 , 
n 2/3 si(x) -> si, u 2/3 si(y) -> si, 
n 2/3 s 2 (x) s 2 , ri 2/3 s 2 (y) ->■ s 2 , 



(8.9) 
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as n — > oo (see Lemma 17. 15[) . Note that r\(x) = r\{y) — r\\ compare (|7.44l) - (|7.45l) with (|2.50l) 
(|2.51[) . Furthermore, we have that 



R~ 1 (y)R(x) =1 + 



x-y 

1/3 



I + O 



(8.10) 



and in view of (|7.49p we see that E n (x) = 0{n^ 6 ), E n {y) = 0(n^ 6 ) and 

E- 1 (y)E n {x) = I + 0{n- 1 ' 3 ), 



ill) 



as n — > oo. The constants implied by the O-symbols are independent of u and v when u and v 
are restricted to compact subsets of the real line. Thus, a combination of (|8.10p - (|8.11l) and the 
fact that N(z) — Nq/z is uniformly bounded with respect to z near the origin and n gives us 



lim E-\y) (l+N(y)-^) R- 1 (y)R(x) (i — N(x) + — ^ E n (x) 
n->oQ \ y J \ x J 

Inserting this into ([57B]) . we then obtain from and (|577| - (|5^| that 

k*(u— v)n 1/3 /c 



= I. 



lim • 

n— f oo 



7 2 /3 



t u v \ 
V cn 2 / 3 ' cn 2 / 3 / 



1 



27ri(u — u) 
This, together with (g^g}, yields 



(-1 1 0) M + 1 {u- 1 r ll r 2 ,s 1 ,S2)Mj r (v\r 1 ,r2,s 1 ,s 2 ) 




1 

\V 



1.12) 



e k'(u-v)n 1/3 /c 

lim — K„ 

n—>oo an, I 



- en 2 / 3 ' cn 2 / 3 



= K 



tacnode 



{u,v;r 1 ,r 2 ,s 1 ,s 2 ), 



U3) 



which is (|2.49p . since c 2 =k*/c. 

The case where x and/or y are negative, or equivalently, u and/or v are negative can be 
proved in a similar manner. We do not give details here. 

This completes the proof of Theorem 12.71 □ 



8.2 Proof of Theorem [2TT41 

In this section we prove Theorem 12.141 The proof will be quite similar to the one in [T2] . 

Following the transformations (|7.3I) . (I7.10[) , (|7.18[) . (|7.57[) and (|7.66l) in the steepest descent 
analysis, we have the following representation for z sufficiently large in the region between the 
contours r r and T;: 

Y(z) = L n U{z)A~ n (z) = L n R{z) (I + N Q /z) p(°°) (z)A' n (z), (8.14) 
where L and A are given in (|7.12l) and (|7.11l) . As in [19], the following lemma is easy to check. 
Lemma 8.1. For the matrix Y± in (12.131) we have 

Yi = L n (Ai + P 1 (oo) + A^o + Ri) L~ n , (8.15) 



GG 



where A\, P^°°^ and Ri are matrices from the expansions as z — > oo, 



K-' n {z)L r ' 



Ai 



i + — + o\^ 



1 



P^(z) = 1+^- 



(oo) 



O 



diag(z™ 1 ,z" 2 ,z-" 1 ,z-" 2 ), (8.16) 

(8.17) 
(8.18) 



We are only interested in the combinations (^1)1,2(^1)2,1 an d (^1)1, 4(^1)4,1 of entries of Y\. 
Since L is a diagonal matrix, the factors L n and L~ n in (|8.15D will not play a role for these 
combinations. Also, since Ai is a diagonal matrix (which is clear from (|8.16[) . since A(z) and L 
are both diagonal), this does not play a role either. Therefore we have for i < j, 



R{z) = I+?±+o(\ 



(>U..iO'i), 



p, 



In what follows we evaluate p.*- 00 -* , Nq and R\ 



(00) 



N Q + Ri 



P^ ] +N + R 1 



(8.19) 



Lemma 8.2. The matrix F 1 (oo) in (p3~TT)) can be written as 



/ 










\ 

















VpI 











V 










/ 



(8.20) 



This result and its proof are exactly the same as in |19j . 
Lemma 8.3. The matrix R\ in (|8.18l) satisfies 

Ri = 0(n~ 2/3 ), as n ->• 00. 
Proof. The asymptotic behavior of Jr in (|7.67l) can be extended to 



Mz) = I + j£\z)n- 1 / 3 + 0(7 



-2/31 



.21) 



3.221 



as n — > 00, uniformly for all z on the circle \z\ — n for a certain matrix-valued function 
JjP(z). An explicit formula for Jr\z) can be obtained in terms of the matrices Mi and M3 in 
(|7.50[) . For us, it will be sufficient to know that this matrix has a Laurent series expansion at 
the origin: 

J Ft\ z ) = ^3,n~^ + Ci n —z + C\ n h Co ,n( z )> 

z 6 z z z 

where the matrix function Co, n (z) is analytic at the origin, and where the matrices Cfc,„, k — 
1,2,3, are depending on n in such a way that the limiting matrices 



lim n 1/3 C k . n =■ C k 



3.23) 



exist for all k = 1,2, 3. 

From (|8.22[) . standard considerations [TTIISI] show that R itself has an expansion of the form 



R(z) =1 + RW(z)n- 1/3 + 0(n- 2 / 3 ), 



3.24) 
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for n oo, where the leading term i?W is given by 

= ! ^ + + if |z|>n-V3, 

t -(-o,n(4 else. 

Taking into account (|8.23l) ~ (|8.25p . we see that the matrix R\ in (18. 18)) must be of order 0(n~ 2 / 3 ), 
and the lemma follows. □ 

Lemma 8.4. We have 

No=K \C d) n_1/3 + °( n_2/3 )' asn->oo, (8.26) 
where K is the constant defined in (|2 . 59[) . and where the 2x2 blocks A, B, C , D are given by 
A = t(bt-b* 2 )( q * a) -f)), (8.27) 

B = iy/t(l-t)(a* 1 -a* 2 )(b* 1 -b* 2 ) (* a) q &Y (8.28) 

C = -iy/t(l - t)(at - a* 2 )(bt - b* 2 ) (* a) q{ ^y (8.29) 

D = (l-t)(at-4)Q a) _ « (<r) ), (8.30) 

with the value a given by (|2.54l) . with q denoting the Hastings- McLeod solution to Painleve II, 
and where the entries denoted with * depend on the Hamiltonian u through the constants c, c in 
(j232jl - ([2^3)1 . 

Proof. Recall that by definition, No is the residue of the matrix N in (|7.54[) at the origin. The 
lemma then follows by computing this residue on account of the formulas in (|2.27l) . (|2.41[) - (I2.43|) . 
(|230]) - (l23Tj) . (j733l and fl7T4"3"]) . taking also into account Remark □ 

Remark 8.5. For facility of comparison, we have stated the above lemma in exactly the same way 
as in PI5] . We note however that the statement can be considerably simplified in the present case, 
since we are now looking exactly at the critical time t = i cr it- This means that the variable t in 
the above formulas can be substituted by (|2.27p . This implies in particular that in (|8.27p - (|8.30p 
we have 

t(b{ - b%) = ^t(l-t)(al-a* 2 )(bl-b*) = (1 - t)(al - a*). 

Further simplification can be obtained by substituting (|2.26p . 

Lemmas 18.21 and 18.41 have exactly the same form as the corresponding results in [TO] . (The 
entries denoted with * in Lemma 18.41 are slightly different when compared to [19] , but this turns 
out to be irrelevant.) Hence Theorem 12.141 (and also the analogous result for the 'diagonal 
recurrence coefficients') follows from the calculations in |19j . 
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